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Abstract. This is the first in a series of papers jZ3l IZ4| on inverse spectral/resonance 
problems for analytic plane domains O. In this paper, we present a rigorous version of 
the Balian-Bloch trace formula jBBll IBB2| . It is an asymptotic formula for the trace 
Tr\Q,Rp{k + irXogk) of the regularized resolvent of the Dirichlet Laplacian of as fc 



oo 



with T > 0. When the support of /5 contains the length of precisely one periodic reflecting 
ray 7, then the asymptotic expansion of TrlQRp{k + ir logfc) is essentially the same as the 
wave trace expansion at 7. The raison d'etre for this approach to wave invariants is that they 
are explicitly computable. Applications of the trace formula will be given in the subsequent 
articles in this series. For instance, in |Z8I IZ4| we will prove that analytic domains with 
one symmetry are determined by their Dirichlet (or Neumann) spectra. Although we only 
present details in dimension 2 for the sake of simplicity, the methods and results extend 
with few modifications to all dimensions. 



1. Introduction 

This paper is the first in a series of articles devoted to inverse spectral and resonance 
problems for analytic domains C with Dirichlet (or Neumann) boundary conditions 
|Z3| IZ4j . As in the earlier articles \Z1\ \Z2\ llSZj , the aim is to recover the domain Q from 
spectral invariants associated to special closed orbits of the billiard flow, specifically the wave 
trace invariants B^j at iterates Y of a single bouncing ball orbit 7. Such wave invariants 
are polynomials in the Taylor coefficients f±\o) of the defining functions f± of at the 

endpoints of 7 and our goal is to recover the coefficients f±\o) from the invariants B^j as 
r varies, and hence to recover the analytic domain. In this series, we will achieve this goal 
when Q satisfies one symmetry condition. Except at the very last inverting step of |Z3| IZ4j . 
we do not use any symmetry or analyticity assumptions on fl and our computations of wave 
invariants are valid on all smooth plane domains. We restrict to plane domains to simplify 
the exposition; the methods extend in a straightforward way to domains in M". 

The path we take towards the wave invariants Bjrj is the one initiated by Balian-Bloch 
in the classic papers [BBl, BB2]. As recalled below, these papers are concerned with the 
asymptotics of a (regularized) resolvent trace rather than the trace of the wave group, and 
are based on the Neumann series representation for the resolvent kernel in terms of the free 
resolvent. Better known to physicists than to mathematicians, the Balian-Bloch papers were 
one of the origins of the Poisson relation for manifolds with boundary (see below and |CV2j ). 
The purpose of this first paper in the series is to give a rigorous version of the Balian-Bloch 
approach. In |Z3j . we use it to calculate wave invariants in terms of Feynman diagrams and 
amplitudes and to reduce the inverse spectral problem to concrete combinatorial problems. 
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In subsequent articles, we will demonstrate the usefulness of this reduction with the fol- 
lowing applications: 

• In jZ8j we prove that simply connected analytic domains with one symmetry (that 
reverses a bouncing ball orbit) are determined among such domains by their Dirichlet 
spectra; 

• In |Z4j we prove that exterior domains with one mirror symmetry are determined by 
their resonance poles. This is a resonance analogue of our inverse spectral results; 

• These results should admit extensions to domains formed by flipping the graph of 
an analytic function around the x-axis, allowing for corners at the x-intercept. It is 
easy to construct examples of this kind. 

The crucial advantage of our present approach is the very explicit nature of the trace 
formula, which allows us to obtain relatively simple formula for lower order terms in wave 
trace expansions. These lower order terms allow us to remove one symmetry from the inverse 
results of ourself and of lantchenko-Sjostrand-Zworski [ ZT\ IZ21 llSZj . The rather concrete 
formulae for wave invariants (see IT^ are valid for all plane domains and the symmetry 
assumption is only used to reduce the amount of data required to determine the domain. 

We should emphasize that we view the Balian-Bloch asymptotic expansions as primarily a 
computational device. Although one could give a new and self-contained proof of the Poisson 
relation for bounded Euclidean domains by this approach, we do not do so here. In particular, 
we use the Melrose-Sjostrand results on propagation of singularities to microlocalize traces to 
orbits. Also, we do not explain how to cancel so-called 'ghost orbits' of non-convex domains, 
i.e. closed billiard trajectories which do not stay within the domain. To explain how the 
Balian-Bloch approach improves on normal forms or microfocal parametrix constructions, 
and to explain its connection to the usual Poisson relation, we now give an informal exposition 
of the Balian-Bloch trace formula. 

1.1. Introduction to wave trace and Balian-Bloch. Let Q denote a bounded C°° plane 
domain, and let 



denote the even part of the wave group of the Dirichlet Laplacian Af^. As recalled in §3, 
the singular points t of the distribution trace TrEQ{t) are contained in the length spectrum 
Lsp{fl) of Q, i.e. the set of lengths t = of closed orbits of the billiard flow of Q 
(i.e. straight-line motion in Q with the Snell law of reflection at the boundary; cf. §2). In 
particular, at lengths of periodic reflecting rays 7 the trace has a complete singularity 
expansion; its coefficients are known as the wave trace invariants of 7. 

The wave trace invariants at a periodic reflecting ray may be obtained from a dual semi- 
classical asymptotic expansions as k —>■ 00 for the regularized trace of the Dirichlet resolvent 



We emphasize that Aq (the usual Dirichlet Laplacian) is a negative operator, so the signs 
of the two terms are opposite. In the classical work of Seeley and others on resolvent traces, 
asymptotics are taken along the vertical line (or a ray of non-zero slope) in the upper half 
plane k + it and traces are polyhomogeneous functions of k. Here, we are taking asjmap- 
totics along horizontal lines (or logarithmic curves) in this half-plane and obtain oscillatory 
asymptotics reflecting the behaviour of closed geodesies. 
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We fix a non-degenerate periodic reflecting ray 7 (cf. §2). For technical convenience we will 
assume that the reflection points of 7 are points of non-zero curvature of dfl. This is a minor 
assumption, but it simplifles one argument f Proposition 17. Ij) and we use it again in |Z3j . We 
briefly explain there how to remove this assumption. We also let p G C^{L^ — e,L^ + e) be 
a cutoff, equal to one on an interval {L^ — e/2, + e/2) which contains no other lengths in 
Lsp(f2) occur in its support. We then deflne the smoothed (and localized) resolvent by 

(1) Rp{k + iT):= p{k - fi){fi + iT)Rn{fi + iT)dfi. 

Jr 

When 7,7~^ are the unique closed orbits of length L^, it follows from the Poisson relation 
for manifolds with boundary (see §3 and |GM1 IPSj : see also Proposition (jSH})) that the 
trace TrlQRp{{k + ir)) of the regularized resolvent on L^(fi) admits a complete asymptotic 
expansion of the form: 

00 

(2) TrlnRp{k + ir) ~ e^*^"")^^ 5Z(57,i + B^-\j)k~^ , k ^ 00 

i=i 

with coefficients B^j, B^-i j determined by the jet of Q at the reflection points of 7, in the 
sense that 

R 

The coefficients B^ j, B^-i j are thus essentially the same as the wave trace coeffflcients at the 
singularity t = L^. Our main goal in this paper is to give a useful algorithm for calculating 
them explicitly in terms of the deflning function of dQ. 

In fact, it is technically more convenient to consider asymptotics of traces along logarithmic 
curves k + ir log k in the upper half plane. We therefore modify the regularization (0) to 

(3) Rp{k + IT log k) := I p{k — p){p + iT\ogk)RQ{p + irlog ix)dp. 

Jr 

In place of © we will get 

00 

(4) TrlnRpik + irlogk) ~ e^^'^^^^fc""^^ + B^~i,j)k-^, k ^ 00. 

i=i 

The additional power law decay k^'^^"' will not cause problems in our study of wave invariants 
at a flxed closed geodesic, because the the errors have the accuracy of k'°°. 

In principle, one could obtain sufficiently explicit formulae for B^j by applying the method 
of stationary phase to a microfocal parametrix at 7 (cf. |GMj |PSj ). or by constructing a 
Birkhoff normal form for A at 7 |Z11 IZ21 ISZ ISZ . However, in practice we have found the 
Balian-Bloch approach more effective. Its starting point is to write the Dirichlet Green's 
kernel Gn{k + iT,x,y) of Rn{k + ir) as a Neumann series (called the multiple reflection 
expansion in |BBlj ) in terms of the free Green's function Go{k + iT,x,y), i.e. the kernel of 
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the free resolvent Rn{k + ir) = — (Aq + {k + ir^) ^ on R^: 

Gn{k + iT,x,y) = Go{k + iT,x,y) + YlZ=i{-'^Y' GM{k + iT,x,y), v^heie 

(5) GM{k + iT,x,y) = J^q^^i d^,^Go{k + ir, x, qi)Go{k + ir, qM,y) 

X Ufjl^d^^Goik + ir, g^+i, qj)ds{qi) ■ ■■ds{qM), 

where ds{q) denotes arclength on dQ and d^^ is the interior unit normal operating in the 
second variable. (The only change in the case of the Neumann Green's kernel is that the 
Mth term has an additional factor of (—1)^; by making this change, all our methods and 
results extend immediately to the Neumann case.) The terms are regularized as in by 
setting 

(6) GM,p{k + ir log k)) = j p{k - ^){^ + ir log ^)Gm{^- + it log ^)d^ 

There exists an explicit formula for Golk + ir) in terms of Hankel functions (§4), from which 
it appears that the traces Trl^^G M,p{k+iT log k) = GM,p{k+iT log k, x, x)dx are (formally) 
oscillatory integrals with phases 

L{x,qi, . . .,qM,x) = \x - qi\ + |gi - ^2! H h \qM - x\, 

equal to the length of the polygon with vertices at points (gi, . . . , g^) £ (dQ)^' and x E Q. 
The smooth critical points correspond to the M-link periodic reflecting rays of the billiard 
flow of fl of length L^, satisfying Snell's law at each vertex (for short, we call such polygons 
Snell polygons). Since the amplitudes and phases only involve the free Green's function, they 
are known explicitly and the coefficients of the stationary phase expansion of Tr l^Rp{{k + 
ir log k)) can be calculated explicitly. 

Thus, our plan for determining Q from wave invariants at a bouncing ball orbit 7 is as 
follows: We chose p to localize at the length of the rth interate of 7. We then localize the 
integrals G {k + ir log k) to small intervals around the endpoints {(0, — L/2), (0, L/2)} and 
parametrize the two components, as above, hj y = f-{x), resp. y = f+{x). (see figure ??) The 
integral over (dQ)^^ is then reduced to an integral over (— e, e)*^ with phase and amplitude 
given by canonical functions of /_(a;), /+(x). Applying the stationary phase method, we 
obtain coefficients which are polynomials in the data f^\o),f^\o). By examining the 
behaviour of the coefficients under iterates Y of the bouncing ball orbit, we try to determine 
all Taylor coefficients f±\o) from these stationary phase coefficients. In jZ3UZ4l[Z5] . we will 
see that this method suceeds at least if the domain has one symmetry (so that /+ = — /-). 

A number of technical problems must be overcome in turning this idea into a proof. The 
problem is that the free Green's kernel (and its normal derivative) only possesses a WKB 
formula away from the diagonals = qi+i, and worse, it is singular along these diago- 
nals. Hence, in the Balian-Bloch approach (unlike the parametrix approach), the following 
problems arise: 

• The individual terms Gm,p niust be regularized, i.e. converted into sums of standard 
oscillatory integrals. Until then, it is not even clear that the trace of each term has 
asymptotic an asymptotic expansion, or that it localizes at critical points; 
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• There are non-smooth critical points (corresponding to Snell polygons in which at 
least one edge collapsed), and their contribution to the stationary phase expansion 
must be determined; 

• Since the multiple reflection expansion (jSJ is an infinite series, it must be explained 
how the full series produces an asymptotic expansion (0) of Trlo^Rpik + irlogA;), 
and why each term in Q depends on only the trace of a finite number of terms of 
(jSJ. In particular, the 'tail' trace must be estimated. 

Let us briefiy outline how we deal with these difficulties. At the same time, we will justify 
the effort involved by explaining how the new results are obtained by this method. 

1.2. Regularizing the terms of the Neumann series. We first explain how we regularize 
the individual terms in the Neumann series. For expository reasons, we suppress the role of 
the interior variable x in the introduction. 

As will be explained in ^ the multiple refiection expansion is derived from an exact 
formula of potential theory (jl^ for the Dirichlet resolvent: 

(7) Rn{k + it) = Ro{k + ir) - Vi{k + it){I + N{k + iT))-hnSt{k + it), 

by expanding the operator (/ + N{k + it))'^ in a geometric series. Here, + it) (resp. 
Sl{k + iT)) is the double (resp. single) layer potential (see (jSH)), S^'^{k + iT) is the transpose, 
and N{k + iT) is the boundary integral operator on L^((9f2) induced by Vi{k + iT) (see 
Also, Ro{k + it) is the free resolvent on M^, and is the restriction to the boundary. The 
existence of the inverse (/ + N{k + iT))~^ is guaranteed by Fredholm theory. 

As will be discussed in detail in ^ (see also H^), the operator N{k + it) is a hybrid 
Fourier integral operator. It has the singularity of a homogeneous pseudo differential operator 
of order —1 on the diagonal (in fact, it is of order —2 in dimension 2, see Proposition 14. 
This is the way it is normally described in potential theory [Tj. However, away from the 
diagonal, it has a WKB approximation which exhibits it as a semi-classical Fourier integral 
operator with phase dgfi^q, q') = \q — q'\ on dQ x dfl, the boundary distance function of Q. 
To make this precise, we introduce a cutoff xi^^^^lQ ~ q'\) to the diagonal, where 6 > 1/2 
and where x ^ C'^(]R) is a cutoff to a neighborhood of 0. We then put 

(8) N{k + iT) = No{k + iT) + Ni{k + iT), with 

r Noik + IT, q, q') = xik'-^\q -q'\)N{k + it, g, q'), 

(9) <^ 

( m{k + IT, q, q') = (1 - x{k^-'\q -q'\)) N{k + it, q, q'). 

As will be shown in Proposition 14. 3[ Ni{{k + iT),q,q')) is a semiclassical Fourier integral 
operator with phase equal to dgn^q, q'). 

Roughly speaking, the boundary distance function generates the billiard map of dVL (see 
^2.21 for the definition). This is literally correct only on convex domains, since dg^{q,q') 
generates both the interior and exterior billiard map, and hence on non-convex domains its 
canonical relation contains 'ghost orbits' (orbits which may exit and re-enter the domain). 
Because we are microlocalizing to one periodic refiecting ray 7, ghost orbits play no essential 
role in this paper and we think of Ni{k + it) as 'quantizing' the billiard map. 
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Now consider the powers A^(A; + ir)^'^ which arise when expanding (J — N{k + ir)) ^ in a 
geometric series. We write 

(10) {No + N,)^= Yl Kii)°N^i2)0---oN^^M). 

a:{l,...,M}^{0,l} 

To regularize N'^ is essentially to remove all of the factors of Nq from each of these terms. 
This is obviously not possible for the term Nq^ but it is possible for the other terms. In 
Proposition 16. 11 we show is that NqNi and NiNq are semiclassical Fourier integral operators 
of the same type as A^^i (and with the same phase), but with an amplitude of one lower degree 
in k. Thus, the term N^^ is of the highest order in the sum. We will see that it is the only 
important term for the inverse spectral problem. In Proposition 17.11 we similarly break up 
the layer potentials in ((7j) and analyse compositions with these. Further we compose with a 
special kind of semiclassical cutoff operator x(x, k~^Drc) to a neighborhood of the orbit on 
both sides of ((Tj). 

1.3. The M-aspect. Since the geometric series expansion of {I + N{k + iT))~^ is very slowly 
converging, we write it as a partial geometric series and a remainder: 

{I + N{k + ^T))-' = J:^;Ui-lrNik + ^Tr^' + 7^Mo, 

(11) 

where TZmo = N{k + ir)*^"+i (/ + N{k + ir))-'^. 

The partial geometric series is regularized by the methods described in the previous section. 
We now explain how to estimate the tail when evaluating wave invariants at a closed orbit 7. 
In our applications, 7 is a bouncing ball orbit, so we will assume it is one in the remainder 
estimate. 

We use elementary inequalities on traces to reduce the estimate of the tail trace 
(12) Trlf,7^Mo {k + IT log k)x{x, k~'D,) 

to an estimate of TrN'^'^Xd-yi^+iT log k)xg^N^'^°* , where xa-y is a semiclassical cutoff operator 
on dfl to the periodic orbit of the billiard map corresponding to 7. As will be verified in 
[|H1 after expanding as in (fTUj) . the operator N^^° Xd-yik + iT log k)xQ^N^'''°* can be regularized 
as above as a sum over a of semiclassical Fourier integral operators whose phases involve 
lengths C of Mq — |cT|-link billiard trajectories, where \a\ represents the number of A^o factors 
in a term. The cutoff operator will force the direction of the first link to point nearly along 
7 and hence will force each link of the critical Mq — |cr|-link paths to have length ~ L^. 
We then use the fact that the spectral parameter varies over a logarithmic curve: due to the 
imaginary part of rlogfc of the phase, each link in the critical path gets weighted by the 
factor e~'^^T^°s'^ = k~'^. On the other hand, each removal of a factor Mq lowers the order in 
k by one. Combinng these two effects, we show in Proposition 18.31 and Lemma (8. II that for 
Mq sufficiently large, the remainder term will be of lower order than any prescribed power 

1.4. Main results. This regularization procedure provides a method of explicitly calculat- 
ing wave invariants associated to an M-link periodic reflecting ray. After regularizing the 
initial terms of the multiple reflection expansion and discarding the remainder, we may apply 
stationary phase to rather simple and canonical oscillatory integrals. 
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Our main result may be summarized in the following theorem. The term 'canonical' means 
'independent of the domain 

Theorem 1.1. Let 'j be a primitive non-degenerate m-link reflecting ray, whose reflection 
points are points of non-zero curvature of dQ, and let p G C^(R) be a cut off satisfying 
supp p n Lsp{Q) = {rL^}. Then there exists an effective algorithm for obtaining the wave 
invariants at 7 consisting of the following steps: 

(i) Each term GM,p{k + irlogk) of ^ defines a kernel of trace class; 

(a) The traces TrlQGM,p{k + irlogk) can be regularized in a canonical way as oscillatory 
integrals with canonical amplitudes and phases; 

(Hi) The stationary phase method applies to these oscillatory integrals. The coefficients B^j 
of a given order j < R are obtained by summing the expansions for TrlQGM,p{k + irlogk) 
for M < Mr. There exists Mq such that the tail trace / fl^) is 0{k~^). 

(iv) The coefficients of the term B^.j are universal polynomials in the 2j-\-2-jet of the defin- 
ing function of dQ at the reflection points. 

(v) If 'J is a bouncing ball orbit, then modulo an error term -R2r(i^"'^^/(0)) depending only 
on the {2j — 2) -jet of f at the endpoints, we obtain the formula 

B,r,_, + B,-r^^_, = a,,.,+/f )(0) + a,,,-f^-'\0) + ^(0) + b,,r,-f^''''\o), 

where the coefficients are polynomials in the matrix elements h^"^ of the inverse of the Hessian 
of the length function L at Y with universal coefficients. 

For instance, if 7 is elliptic and invariant under an isometric involution a of f2, then 
(13) 

1 

^ ^ COS Ol I ^ 

q=l 

Here, e^^'^*" are the eigenvalues of the Poincare map of 7. To our knowledge, no such 
explicit formula for the wave invariants has been available before. 

The most important aspect of Theorem (ll.lj) is the algorithm for computing the coefficients 
of B^j in terms of the defining functions of (9f2, as illustrated in (v) and by the explicit 
formula ()13|) . That is truly what distinguishes the result from previous approaches to 
wave invariants and the Poisson relation. For instance, recovering a Z2-invariant dVL is 
tantamount to sifting the Taylor coefficients of / at from the expression ()13|). and we see 
that this amounts to a fine study of sums of powers in the matrix coefficients k^^ of the inverse 
Hessian of £. We needed the complete symbol of the wave trace to obtain this result, and no 
such details have previously been obtained by any other method. The gain in effectiveness 
is borne out in the subsequent articles |Z3t IZ4j . Although we will not do so here, the proof 
extends with no essential changes bounded smooth domains in R". 

The author would like to thank A. Hassell for collaboration on a related project [HZ] which 
clarified many aspects of the boundary integral operators. He would also like to thank J. 
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to thank Y. Colin de Verdiere for many comments, corrections and constructive criticisms on 
earlier versions of this article, both in its written form and in various verbal presentations. 
His encouragement and insight during the course of this project was much appreciated. 



Contents 



1. Introduction! 1 



1.1. Introduction to wave trace and Balian-BlochI 2 



1.2. R,egularizing the terms of the Neumann seriesi 5 

3 



1.3. The M-aspectl 6 



1.4. Main resultsi 6 



2. Billiards and the length functional! 
2.1. Length functional! 9 



2.2. Billiard flow and length spectrum I 10 

3. Wave trace and resolvent trace asvmptoticsl 11 



3.1. Resolvent and Wave groupl 11 



3.2. Wave trace and resolvent trace asvmptotics! 13 



3.3. Microlocal ciitof!l 14 



4. Multiple reflection expansion of the resolventi 16 



4.1. Multiple-reflection expansion: Dirichlet boundarv conditiona 1^ 



4.2. The operators and NA 20 



4.3. Laver potentials as semi-classical Fourier integral operatord 22 

4.4. Integral formulad 24 



4.5. The case of the unit disd 25 



5. Microlocalizing the trace! 26 



6. Regularizing the boundarv integrals I 27 



6.1. The compositions A^n o A^il 28 



6.2. Proof of (A)\ 32 



6.3. The term o y^J 33 



7. Regularizing the interior integral I 34 

8. Tail estimate 40 

9. Compl etion of the proof of Theorem O 44 
R.eferences! 47 



2. Billiards and the length functional 

We collect here some notation and background results on plane billiards which we will 
need below, mainly following the reference Kozlov-Trechev |KTj . 

Let Q denote a simply connected analytic plane domain with smooth boundary dQ of 
length 271. We denote by T = M\27rZ the unit circle and parametrize the boundary counter- 
clockwise by arc-length starting at some point go ^ dQ: 



(14) q:T^dQcR', = |g(0)| = 1, g(0) = go- 
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We similarly identify the m-fold Cartesian product (dQ)"^ of dQ by T™', and denote a point 
of the latter by (0i, . . . , 0m)- 

By an m-link periodic reflecting ray of Q we mean a billiard trajectory 7 which intersects 
dfl transversally at m points g(0i), . . . , of intersection, and reflects off dfl at each 

point according to Snell's law 



- 9(0i) - 



Here, Ugi^^) is the inward unit normal to dQ at g(0). We refer to the segments q{(pj) — 

as the links of the trajectory. An m-link periodic reflecting ray is thus the same as an m-link 

polygon in which the Snell law holds at each vertex. Since they will come up often, we make: 

Definition 2.1. By -P(0i,. ..,</,„) we denote the polygon with consecutive vertices at the points 
(g(0i), . . . , q{4>m)) £ (dQ)"'. The polygon is called: 

• Non-singular if (f)j ^ for all j 

• Snell if P{(j)i,...,(j,^) is non-singular and if / f73j) holds for each pair of consecutive links 

{^('/'i) - - q{<f>j)}; 

• Singular Snell if P{(f,i,...,(i,^) has fewer than n distinct vertices, hut each non-singular 
pair of consecutive links satisfies Snell's law. 

We will denote the acute angle between the link — and the inward unit normal 

^gi4>j+i} by — <l{(pj),J^q((f,j)) and that between — q{(pj) and the inward unit 

normal at g(0j) by Z(g(0j+i) - q{(t>j), z/g(<^.), i.e. we put 

na\ Qi^j+i) - Qi^j) /II, N / , N 

(16) \q{<\>,^.)-qm ■ "^^^^^ = - ^('^^■)' -^(^.)- 

We also use the notation Z(g(0j+r— g(0j), i^q(^j)) for the angle between the link q{<t>j+r)~q{.4'j) 
and the unit inward normal at q{4>j). 

The function Z(g(0j+i) — q{(pj), ^q{<f>j)) is well-defined on T"^ minus the diagonals Aj^+i = 
{(pj = 0j+i}. It has a continuous extension across the diagonals according to the following 

Proposition 2.2. cosZ(g(0) - g(0')), z/g(^)) = -iK(0)|0' - 0| + O(|0' - 0^). 

Proof (cf. IE] and jZl] ) We have: 

1 

2' 

Now divide by 1^(0) — q'(0')l- 

□ 



(g(0')) - g(0)) ■ ^,(0) = -7t(</' - ^O'/^l^) + O((0 - 0'^^^ 



2.1. Length functional. We first define a length functional on T by: 



(17) L(0i,...,0A/) = |g(0i)-g(02)| + ...+ |g(0^,_l)-g(0Af)|. 
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It is clear that L is a smooth function away from the 'large diagonals' ^jj+i '■= {4>j = 
where it has |x| singularities . We have: 

^\<li(Pj) - = -sinZ(g(0j) - z/g(0^)), 

^\<li(pj) - = sinZ(g(0j+i) - q{(pj), i^g(^^_^,)) 

=^ Wj^ = sinZ(g(0j+i) - q{(i)j),Vq(^^+,)) - sinZ(g(0j) - 

The condition that -^L = is thus that the 2-link defined by the triplet q{4>j), Qi+i) 

is Snell at A smooth critical point of L on is thus the same as an M-link Snell 
polygon. 

We will also be concerned with the length functional L : Q x T*^ defined by: 

(19) 

L(x,0i, . . . ,0Af) = \x- g(0i)| + |g(0i) - g(02)| H h |g(0j\/-i - g(0Af)| + Mm) - x\, 

which is smooth away from the diagonals x = g(0i),x = q{<t>M) together with Ajj+i. Its 
gradient in the x-variable is given by 

(20) V.L = ^^^4^ + ~ '^'^''^ 



|x-g(0i)| |x-g(0M)r 

so that a smooth critical point x of L(x, 0i, . . . , 0m) corresponds to triple (0m, a;, 0i) whose 
2-link is straight at x. We sum up in the following well-known proposition, due to Poincare. 
For background, see |KTj . 

Proposition 2.3. A smooth critical point (a;, 0i, . . . , 0m) of L on Q x T*^ corresponds to 
an M-link Snell polygon with vertices (x, 0i, . . . , 0m)- 

We will also need the formula for the interior normal derivative = i'q((f)^^-^) ■ Vy along 
the boundary of the link-lengths : 

(21) - y\y=<i{^,+i) = igS'jlgi^^t!)! ' = cos Z(g(0,-+i) - g(0,-), z/g(<^,.)) 

2.2. Billiard flow and length spectrum. The (geometer's) billiard flow $* of is the 
flow on T*Q which is partially defined by Euclidean motion in the interior and Snell's law of 
reflection at the boundary. We refer to the billiard orbits as trajectories or rays, and when 
they have only transversal intersections with the boundary we refer to them as transversally 
reflecting rays. The straight line segments between intersection points are called links. 

At (co-) vectors tangential to dQ this law does not uniquely define the flow unless Q is 
convex, in which case tangentially intersecting rays can only travel along the boundary. In 
the non-convex case, there exist rays which intersect dQ tangentially and at such points the 
geometer's billiard flow is not uniquely defined. The propagation of singularities theorem for 
domains with boundary (cf. |AM|, Im51 IPSj ) largely resolves this ambiguity, and completely 
resolves it for analytic domains. It defines the billiard flow as the broken bicharacteristic 
flow of the wave operator, i.e. the trajectories along which singularities of solutions of the 
wave equation move. Roughly speaking, the trajectories are transversally reflecting rays and 
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limits of such rays with many small links. The limit rays intersect the boundary tangentially 
and then glide for some time along the boundary and then re-enter the domain. Since small 
links can only occur in the convex part of the boundary, the entrance and exit points to the 
boundary of a plane domain occur at its inflection points. There exists unique continuation 
of geodesies unless dQ has infinite order contact with a tangent line, and of course this cannot 
occur for analytic domains. For background and further discussion we refer to |PSt IGM| IM] . 
Pictures of gliding rays may be found [GM\ IM] . 

In the Poisson relation for the wave equation, it is the analysts' billiard flow (propagation 
of singularities) which is relevant, and henceforth we assume the billiard flow defined as 
the broken bicharacteristic flow for the wave equation. We then define the length spectrum 
Lsp{Q) to be the set of lengths of periodic orbits of the billiard flow ( |PSj . Definition (1.2.9)). 

By the billiard map jS of fl we mean the map induced by $* on B*dfl: if (g, r]) G B*{dQ), 
we may add a multiple of the unit normal to obtain an inward pointing unit vector v at q. 
We then follow the billiard trajectory of v until it hits the boundary, and then define (3{q, rf) 
to be its tangential projection. 

In the case of strictly convex domains, periodic orbits are either periodic reflecting rays 7 
or closed geodesies on dVL { |PSj . Ch. 7). By periodic n-link reflecting ray, we mean a periodic 
orbit of the billiard flow $* on T*Vl whose projection to Vl has only transversal intersections 
with dVL. That is, 7 is a Snell polygon with n sides. (Here, and henceforth, we often do not 
distinguish notationally between an orbit of $* and its projection to Vt.) 

2.2.1. Poincare map and Hessian of the length functional. The linear Poincare map of 
7 is the derivative at 7(0) of the first return map to a transversal to at 7(0). By a 
non-degenerate periodic reflecting ray 7 we mean one whose linear Poincare map P^ has no 
eigenvalue equal to one. For the definitions and background, we refer to |PSj |KTj . 

There is an important relation between the spectrum of the Poincare map P^ of a periodic 
n-link reflecting ray and the Hessian Hn of the length functional at the corresponding critical 
point of L : T" R. For the following, see |KTj (Theorem 3). 

Proposition 2.4. We have: 

det(/ - P^) = - det{Hn) ■ (61 ■ ■ ■ bn)-\ 

where h, = ^'''^^t^'l,^'^'^^'^' • 

3. Wave trace and resolvent trace asymptotics 

The spectral invariants we will use in determining Q are essentially the wave trace in- 
variants associated to bouncing ball orbits. As will be recalled below, these invariants are 
coefficients of the asymptotic expansion of the trace Trl^E^{t) of the Dirichlet wave group 
around in singularities at lengths of periodic billiard trajectories. Dual to the wave trace sin- 
gularity expansion, at least formally, is the asymptotics as — >■ cx) of the trace Trl^RQlk+ir) 
of the Dirichlet resolvent. Further, we will regularize the trace and relate the resolvent trace 
coefficients at a periodic reflecting ray to the corresponding wave trace invariants. 

3.1. Resolvent and Wave group. By the Dirichlet Laplacian Af^ we mean the Laplacian 
^ = + -§-2 with domain {u E Hq^VI) : Au G L^}; thus, in our notation, is a negative 
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operator. We denote by x, y) = cost^/—AQ{x, y) the fundamental solution of the mixed 
wave equation with Dirichlet boundary conditions: 



(22) 



{ 9^ = AE onRxnxn 

En{0,x,y)=6{x-y) ^(0,x,i/)=0 
^ En{t, x,y) = {t, x,y) eRx dVl x il. 



The resolvent of the Laplacian A^^ on Q with Dirichlet boundary conditions is the operator 
on defined by 

Rn{k + it) = -{An + {k + itY)-\ r > 0. 

The resolvent kernel, which we refer to as the Dirichlet Green's function Gn{k + iT,x,y) of 
f2 C M^, is by definition the solution of the boundary problem: 



(23) 



(A^ + {k + iTf)Gn{k + ir, x, y) = -6{x-y), {x, yen) 
Gfiik + IT, x,y) = 0, X E dQ. 



To clarify our sign conventions, let us specify them in the case of the free Laplacian Aq 
on M^. Our A is negative, so the symbol of the free resolvent is (|^p — {k + ir)"^)^^ . The free 
Green's function has the asymptotic behaviour 

^i{k+iT)\x-y\ 

Go{k + tr, X, y) ~ p ^ ^^^|^ _ ' ([(^ + ' ^0 ^ 

hence is oscillatory in k and has the exponentially decay e""^'^"^" when r > 0. We will later 
place the spectral parameter on the logarithmic curve k + ir log k, where it has a power law 
decay in k as well as oscillatory behaviour. For further discussion of the signs, see [T] (p. 
142). 

The resolvent may be expresed in terms of the (even) wave operator as 



(24) 



Rn{k + it) = —— / e^('=+*^)*En(t)t^t, (r > 0) 



which holds because 

f oo 



1 



X^-{k + zr)2 k + iT 



poo 

/ e^(^+''")* cos Xtdt, (VA eR, re 
Jo 



Given p G C^(M+), we have defined the smoothed resolvent Rp{k + ir) in ((T)). By 
we can rewrite it in terms of the wave kernel as: 

Rp{k + tT) = f^pik- fi)e'^^'+'^^'Enit)dtdfi 

(25) = /o°°p(t)e*('=+-)*i?n(t)rft 

= p{k + iT + y/A^) + p{k + iT- ^/A^) 
This is essentially the smoothing used in the study of wave invariants in jD 
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3.2. Wave trace and resolvent trace asymptotics. We now recall the classical results 
about the asymptotics of TrlnRp{k + it) (cf. |GMj |PSj ). From the last formula in we 
note that 

oo 

(26) TrlnRpik + ir) = ^[p(A; - A^- - ir) + p{k + A^- - tr)]. 
Note that p is an entire function since p G C^(]R). We also note that 

oo 

(27) TrluRp{k + iT) = J2pik-Xj-iT) + 0{k'^), (fc^oo,r>0) 

i=i 

since Xj > for all j. This is essentially the same expression studied in |DG| IGMj . 

Dual to the sum is the trace of the even part of the Dirichlet wave group, i.e. the distri- 
bution in t defined by 

„ oo 

TrlnEQ{t) := / E^{t,x,x)dx = y^costAj 
Jn 

where 

A0J = Aj0j, 0j|an = O, {(j)i,(f)j) = 6ij, En{t,x,y) = ^costXj(j)j{x)(t)j{y). 

j 

The singular support of the wave trace is contained in the set Lsp{Q) of lengths of gen- 
eralized broken geodesies ( \AM\ lUM] . Theorem; and |PSj ) : More precisely, for any bounded 
smooth domain, we have 

singsuppTrluEnit) C Lsp{Q). 

If Q belongs to a certain residual set 71, then 

singsuppTrlnEnit) = Lsp{Q). 

When is the length of a non-degenerate periodic reflecting ray 7, and when is 
not the length of any other generalized periodic orbit, then TrlQE^{t) is a Lagrangean 
distribution in the interval {L^ — e,L^ + e) for sufficiently small e, hence Trl^Ep{k + ir) has 
a complete asymptotic expansion in powers of k~^. Let us recall the precise statement (see 
[(ilMj . Theorem 1, and also page 228; see also jPSj Theorem 6.3.1). 

Let ^ be a non- degenerate billiard trajectory whose length is isolated and of multiplicity 
one in Lsp{Q). Then fort near L^, the trace of the wave group has the singularity expansion 

00 

TrlnEnit) ~ a^{t - + iO)'^ + a^o log(t - + iO) + ^ a^fc(t - + iO)'' \og{t - + iO) 

k=l 

where the coefficients a^^ (the wave trace invariants) are calculated by the stationary phase 
method from the Lagrangean parametrix E. 

Recall that E{t) is a microlocal parametrix in that it approximates E^it) modulo regular 
kernels in a sufficiently small conic neigbhorhood of R"'"7. 
We will need the following equivalent statement: 
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Corollary 3.1. Assume that 7 zs a non-degenerate periodic reflecting ray, and let p G 
C^{Lj — e,L^ + e), equal to one on {L^ — e/2,L^ + e/2) and with no other lengths in its 
support. ThenTrl^Rp{k + iT) admits a complete asymptotic expansion of the form (0). The 
coefficients B^.j are canonically related to the wave invariants a^.j . 

Proof. By (j25|) . we have Rp = Ep where 

(28) Ep{k + iT):= [ e^'^-^^'p{t)En{t)dt. 

Jr 

The corollary thus follows immediately from the Poisson relation. 

□ 

As mentioned in the introduction, we actually use a variant of this result: 

Corollary 3.2. Under the same assumptions, TrlQRp{k + irlogk) admits a complete 
asymptotic expansion of the form ^ with the same coefficients B^j . 

Proof. We use (PHj) but with r log k in place of r. We then subsitute the microfocal parametrix 
E(t) and calculate 

Jr 

asymptotically by the stationary phase method. Since supp p is contained in R_|_, the factor 
g-riogfci j-Qg^y absorbed into the amplitude and decreases its order. The result follows 
exactly as in Corollary ()3.1|) by the stationary phase method. 

□ 

3.3. Microlocal cutoff. In this section, we explain that the regularized wave trace ex- 
pansion TrlQRp{k) at a periodic reflecting ray 7 can be microlocalized to 7, i.e. equals 
TrlRp{k)x{k), where x(^) is a (specially adapted) cutoff to 7. In principle this should be 
obvious, but we include some details since we were unable to find a suitable reference. For 
simplicity we assume that 7 is a bouncing ball orbit. 
The cutoff consist of there terms, 

Xik) = x+{k) + Xoik) + x-{k) 

corresponding to the top, middle and bottom of the orbit. Let f/ = f/+ U [/q U f/_ be a small 
strip around vr(7) in R^, with U± be a small neighborhood oi U Ci dQ± (the top/bottom 
boundary component). Further, let {6,r) — > q{6) + rz/g(6i) denote Fermi normal coordinates 
along U+ and we denote the dual symplectic coordinates by Pr,Pe- We use the same notation 
for f/_, anticipating that no confusion will arise. We then define semiclassical pseuoddiffer- 
ential cutoff operators of the form: 

• X±{''^y^yk~^Dg) on U±, with X+{^^(pyPe) supported in f/+ x {\pg\ < e} and with 
X±{r,e,pe) = 1 for all < r < e,0 e (-e,e), \p^\ < e/2. 

• xo{^: k~^Dx) is properly supported in Uq x Uq and Xo{x,^) is supported in a small 
neighborhood of Uq x {(0, 10}. 

One can apply x{k) to functions supported in f/ fi fi. In particular, the cutoff resolvent is 
defined near the boundary by 

X±{k)RnAk + ^r){r, 9, r', 9') = [ [ e''^'-'"^Po Rn{k + it, r, 9"- r', 9')x±{r, 9,pe)dped9". 

Jr Jr 
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Naturally, the introduction of the cutoff does not change the trace modulo negligeable 
terms. For the sake of completeness, we sketch the proof that the full resolvent trace is 
unchanged. 

Lemma 3.3. TrlnRn,p{k + ir log k) - Trlnx{k)RnAk + ir log k) = 0{k-°°). 
Proof. It is known jXM] (see also [HMj l that 

TrlnRn,p{k + it) ~ TrlnRn,p{k + ir) 

where Rn,p{k + ir) is a microlocal parametrix for the semiclassical Dirichlet resolvent in 
neighborhood of 7. To be precise, Andersson-Melrose |AMj proved that there exists a mi- 
crolocal parametrix for E£){t), the even Dirichlet wave group, near any transversal reflecting 
ray. The Fourier-Laplace transform of this parametrix is a semiclassical resolvent parametrix. 
The parametrix is defined throughout U. 

We verify this using ()25p. First, we may write 

POO 

(29) RpDik)xik)= / pit)EDit)xir,y,\Dt\-'Dy)e'^''+^^^'dt. 

Jo 

Calculating the expansion in the statement of the Lemma is the same as computing the wave 
front set of the trace of ()29j) near rL^. Thus, the statement is equivalent to saying that 

(30) WF[TrED{t){I - x(r, y, \Dt\-^Dy)] n (rL^ - e, rL^ + e) = 0. 
To prove this, we recall that WF{Er)(t,x,y)) is the space-time graph 

(31) ^ = {(^,r,x,e,x,r)e^*(Rx^]^x^]^): r = -|e|, G\x,o = ix',a}, 

of the generalized billiard flow. Now in Fermi coordinates {r, pdr,y,ridy),W~^'j is ray in the 
direction of dr. Since the space-time graph of the cotangent bundle along 7 may be described 
in normal coordinates as a neighborhood of 

{(t,r,t,r,0,0)}, 

a conic neighborhood may be described in these coordinates by \y\ < e, \r]/T\ < e. This is 
precisely the set to which x-yi'^^^y^ \Dt\~^Dy) microlocalizes. Emptyness of the WF in (jHUj) 
follows from the calculus of wave front sets, which implies that only diagonal points in F 
contribute, i.e. periodic orbits of the billiard flow, and from our assumption that Y is the 
only orbit with period in the given set. 

□ 

The only potentially confusing issue is in the choice of cutoff operator, which must be 
rather special since it operators on a manifold with boundary. Let us verify in another way 
that this kind of cutoff operator acts as a microlocal cutoff to 7. 

According to [HMllMlj . we can calculate the trace using a microlocal parametrix 



Eit,x,y)= [ e^^(*'^'^'«)a(t,x,i/,Ocie 

JlR2 
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for the Dirichlet wave kernel. We take its Fourier-Laplace transform to get a semiclassical 
microlocal parametrix 

RnAk + ir. x,y) = / p(t)e^('=+^-)*e^^(*'^'^-^)a(t, x, y, 

JM2 Jo 

We now change variables — > A;.^ to obtain 

RnAk + ^r, x,y) = k^ / p(t)e'^<^(*'^'^'^)e-^*a(t, x, y, k^d^, 
where the phase is $ = t + (j){t, x,y,C,). We now apply x(/c) to get 

Jm? Jo 

where 

d^^) = x{x, d^(f))A{k, X, d^cfy 

with A a symbol of order 0. The assumption on the cutoff implies that x = 1 near 7. Since 
the trace is computed by applying stationary phase to the trace of this oscillatory integral, 
it is unchanged modulo rapidly decaying errors by the cutoff. 

4. Multiple reflection expansion of the resolvent 

The purpose of this section is to review the 'multiple-reflection' expansion of the Dirichlet 
Green's function of a bounded plane domain. This is the term in |BB11 IBB2j for the Neu- 
mann series expression for the Dirichlet Green's function in terms of double layer potentials. 
The same method also works for Neumann boundary conditions, but for simplicity we only 
explicitly treat the Dirichlet case. We refer to jTj, Chapter 5 for background in potential 
theory. 

The method of layer potentials (^T^ II, §7. 11) seeks to solve (j^Hj) in terms of the 'layer 
potentials' Goik + iT,x,q),duyGo{k + iT,x,q) G V'{Q x dQ), where u is the interior unit 
normal to Q, where d,y = u -V, and where Go{k + iT,x,y) is the 'free' Green's function of 
M^, i.e. of the kernel of the free resolvent — (Aq + {k + ir)"^)^^ of the Laplacian Aq on M^. 
The free Green's function in dimension two is given by: 

Goik + IT, X, y) = H^'\ik + tT)\x -y\) = J^, e^(-s/.?)(|^|2 _ + ,^)2)-irf^. 

Here, H^'^^ (z) is the Hankel function of index 0. In general, the Hankel function of index u 
has the integral representations (jTj, Chapter 3, §6) 

(32) 

_ / 2 xi/2 e'("-""/2-"/4) TOO i/2fi s_V-l/2^„ 

~ ^nz) r{u+l/2) Jo ^ ^ 2iz) 

From the first, resp. second, representation we derive the asymptotics: 

' — 2^ ln(|A; + ir| |x — ?/|) as |k + ir| |x — y| ^ 0, if z/ = 



(33) Hl'\{k + zT)\x-y\)r^ < 



— ijirr^^ 1)'' as |k + ir||x-y| ^0, if z/ > 

^i((k+ir)\x-y\-U7r/2~7./4 1 ^ j^l | |x - y| ^ OO. 

(\k+iT\\x—y\)'-/ I ' ' ' ' 
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Here it is assumed that r > 0. 

The single layer, respectively double layer, potentials are the operators 
(34) 

Si{k + iT)f{x) = Jq^ Go{k + IT, X, q)f{q)ds{q), Vi{k + iT)f{x) = jg^ ^Go{k + ir, x, q)f{q)ds{q), 

where ds{q) is the arc-length measure on dQ. They induce boundary operators 
(^) S{k + iT)f{q) = Go{k + IT, g, q')f{q')ds{q'), 

(35) 

(zz) N{k + zr)/(g) = 2 ^G,{k + ^r, g, q')f{q')ds{q') 

which map H'^{dQ) — > H'^^^{dQ). Furthemore one has (cf. jljll. Proposition 11.5) 

(i) S{k + ir), N{k + ir) G ^-\dn), 

(36) 

(ii) (I + N{k + ir)) : H'^^d^l) ^ H'^i^d^l) is an isomorphism. 

To understand the role of these operators, it helps to recall that the Poisson integral 
operator 

(37) Plnik+tr) : H'idn) ^ H'+'/\n), PI^{k+iT)u{x) = I d,Gn{k+iT,x,q)u{q)ds{q) 

J an 

may be expressed in the form 

(38) PIn{k + ir) = 2Vi{k + iT){I + N{k + iT))-\ 

We refer to jT|l (see Chapter 5, Proposition 1.7; see also |LMj ) for background on layer 
potentials. 

We include some other facts about the kernel N{k + ir) which will be used to estimate 
remainders in the multiple reflection expansion (see Proposition ()4.2|1 and Lemma ()8.2|1 . 
Further and more precise estimates will be given in jZ5j . 

Proposition 4.1. Suppose that dVl is Then: 

(i) Nik + IT, g(0i), g(02)) e Ci-(T X T). 

(ii) N{k + ir) is a Hilbert- Schmidt operator on L'^{dQ), with \ \N{k + iT)\\Hs < C|A;|^/^. 

(iii) N{k + ir) G \l'"^(9i7), hence it is a trace-class operator on L'^{dQ). 

Proof, (i) By definition, 

N{k + iT,q{(f)i),q{(f)2)) = d^^Go{k + iT,q{(f)i),y)\y=g(^^^) 

= -{k + iT)H[^\\k + iT\\q{(j)i) - g(02)|) cosZ(g(02) - q{(j)i) , Uqi^^)) . 
Now at r ~ 0, we have (see also [XS], 9.1) and |EPj), 

(39) H[^\r) = — + 0(r(l + \ogr)). 

Tir 

We correspondingly define 

N{k + ir) = Nsing + Nregik + ir), with 
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In fact, the vanishing of the numerator to order 2 imphes that Nsmg is a C°° kernel. Hence 
the smoothness of N{k + ir, q{4>2)) equals that of Nregik + ir, g(0i), q{(p2))- It follows 

that the kernel has the regularity of xlogx, so N{k + iT, g(0i), q{(p2)) is Lipschitz continuous 
of any exponent a < 1. □ 
(ii) Obviously the kernel is Hilbert-Schmidt. We need to estimate 



(40) / \N{k + tT,qi,q2)\'ds{qi)ds{q2). 

For the sake of brevity, we only sketch the norm estimate, referring to )Z5j for further 
details. We break up the domain of integration into the three regions 

\q - q'\ <\k + ir|"\ \k + ir|"^ < |g - q'\ < \k + ir]"^/^, \q - q'\ > \k + ir|"^/^. 

In the first (near diagonal) region, the kernel is bounded above by a constant independent 
oi k + ir. Indeed, both the singularity and the factor of {k + ir) cancel in the most singular 
term, as can be seen from ()39|) or (jHHj) . The smoother terms also cancel the factor of (k + ir) 
in region (i). 

In region (ii), the Hankel factor has a uniform upper bound. The cosine factor puts in 
\q — q'\. Hence, in any region \k + ir|^^ ^ \q — q'\ ^ \k + iT\~^ the integral is dominated by 

"|fc+jr|-'' 

\k + iT\^ I X^dx = 0(\k + tT\'^-^'). 





So it is uniformly bounded if r > 2/3. 

In region (iii), we can (and will) use the WKB expansion of the Hankel function. The 
square of the Hankel function contributes the factor \{k + zr)|~^, cancelling one power in 
(jiUj) . Taking the square root of the result gives the stated estimate. □ 
(iii) It is standard that N{k+iT) G '^^^{dfl) PP. The smoothness of A^smg noted in (i) shows 
that the symbol of of order —1 vanishes; hence, it is actually of order —2. Therefore, it is 
of trace class. 

An alternative proof is by the Hille-Tamarkin theorem (cf. IGDj) theorem, which states 
that a Hilbert-Schmidt kernel K{x, y) on an interval (or circle) is trace class if it is Lipa in 
one of its variables with a > 1/2. □ 

This completes the proof of the Proposition. 

□ 

In |Z5j . we prove that the operator norm of N{k + ir) on L^(f2) is uniformly bounded in 
k for r > 0. On the other hand, it appears that N{k + ir) is not a contraction for any fixed 
r > 0. 



4.1. Multiple-reflection expansion: Dirichlet boundary conditions. To solve 
one puts r(/c + ir, x, y) = Gn{k + ir, x, y) — G^ik + ir, x, y) and solves the boundary value 
problem 

r (A^. - (A: + ir)2)r(A; + zr,x,y) = 0, x,?/Gfi 

(41) I 

y rnT{k + ir, q, y) = -r^Goik + ir, q,y), qe dVL. 

Here, 

r^u = u\on 
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is the restriction operator acting in the first variable. 

We try to solve (PT|) with r{k + ir, y) of the form: r(A; + ir, x, y) = Vi{k + ir)/i. If we 
denote by 

U{q)= hm /(x)GC(n), 

then one has (cf. jTjiI, Proposition 11.1) 

[Vi{k + ir)f] + {x) = l/2/(x) + l/2N{k + iT)f{x). 
Hence the equation for fi{k + ir, q, y) is given by: 

(42) {I + N{k + iT))^ji{k + iT,q,y) = -r^GQ{k + iT,q,y), {q,y) e dQ x Q. 
By (jnni) there exists a unique solution, 

(43) fi{k + iT,q,y) = -{I + N{k + iT))-'rnGo{k + tT,q,y), {q,y) e dQ x Q, 
where by (/ + N{k + ir))-^ G ^^^dQ). It follows that 

R^{k + iT) = Ro{k + iT) -V£{k + iT){I + N{k + iT))-^rnRoik + iT) 

(44) 

= Ro{k + zt) + Em=o(-2)'''^^^(^ + ^r)N{k + irY^mRoik + it). 

More precisely, the kernels are related by the Neumann series 

(45) 

Gn{k + ir,x,y) = Y.Z=^{-'^)^GM{k + ir,x,y), where, for M > 1, 

GM{k + ir, x, y) = J^g^^„ i;Go{k + ir, x, qi)UfL^^-§;^Go{k + ir, qj, qj+i)Go{k + ir, qM, y)UfL^ds{qj) . 

Here, ^ is short for This expansion of Gn{k+iT, x, y) of Rn{k + iT) is referred to in |BBlj 

as the multiple-reflection expansion. We regularize Gm to Gm,p as in (jU)). By substituting in 
the multiple reflection expansion (j45p . we obtain a multiple reflection expansion for Rp. 

In the following proposition, we assume the boundary is smooth since we are dealing 
with analytic boundaries. Some additional work is needed if the boundary is only assumed 
piecewise smooth, as in |Z4j . 

Proposition 4.2. Suppose that dQ is . Then, for r > 0, GM,p{k + iT,x,y) defines for 
each M a trace class operator on L'^{Q). 

Proof. First, the M = term p{k — + iT)RQ{fi + iT)lQdfi is easily seen to be trace 

class since it is the restriction of a smoothing operator to Q. 

Next, we observe that for M > 1, Guik + iT,x,y) is the Schwartz kernel of Vi{k + 
iT)N{k + iT)^~^rQRo{k + ir)!^. Let us consider each factor. Using (jHHj) it follows that 

Vi{k + iT) : H'{dQ) H'+^^^{Q) continuously. 

Further (see T , Chapter 4, Proposition 4.5), the restriction operator satisfies tq : H^{Q) — > 
H'-^/\dQ), for s > 0. (See also [LMl)- And InRoik + ir)ln : H'{Q) H'+\Q). 
It follows that for M > 2, 

Vi{k + iT)N{k + iT)^^-^rnRo{k + tT)ln : H%Q) H'+^^^'^-^'>+\Q). 

Hence, this operator is certainly of trace class on L'^{dQ) if M > 2. 




20 STEVE ZELDITCH 

The case M = 1 is different from the others, so it seems worth considering it separately. 
The operator Vi{fi)rQRo{^)l^ has Schwartz kernel 

p{k - + iT)du^Go{p, X, q{(f)))Go{^i, q{(f)), y)d(f). 

an 

The double layer potential has kernel 

(46) d^^Go{fi,x,q{(f))) = -{iJ, + iT)H[^\{iJ + iT)\q{(f)i) - x\) cos Z{q{(f)) - x,iyg(^^)). 
and as for N{k + ir) above, we write 

duyGoik + ir, x, q{4>)) = Kreg{k + ir, x, g(0)) + Ksing{k + ir, x, g(0)) 

where 

Ksing{k + IT, X, q{(j))) = 2^ [^^Wl iHwi ■ ''"W- 

It is clear that 

/ KregijJ' + ir, X, g(0))Go(/i + ir, q{(j)),y)d(j) 
Jan 

is a composition of a Hilbert-Schmidt kernel from L^(fi) — and one from L'^{dfl) 
L'^{Q). Hence, this term is of trace class. The Kging term has kernel 

1 X — q{(f)) 



an 



|x-g(0)| |x-g(0)| 



J^qWi / p{k - iJ,){iJ. + iT)Go{iJ, + iT,q{(f)),y)dn}d(l) 



is of trace class. However, the bracketed operator is a smoothing operator from L'^{Q) — > 
L'^{dQ). Hence, the composition is of trace class. 

□ 

4.2. The operators A^^o and Ni. We now describe Ni[k + iT) of (jH)) - © as a semiclassical 
Fourier integral operator. For related results, see |HZj . We denote by x a smooth bump 
function which is supported on [—1,1] and equals 1 on [—1/2,1/2]. We also now put in 
h = k-^+^. 

We denote by 5'^(T'") the class of symbols a{k, 0i, . . . , satisfying: 
(47) \ik-'D^raik,<P)\ < (|A;| > 1). 

Proposition 4.3. Ni{k+iT) is a semiclassical Fourier integral operator of order associated 
to the billiard map. More precisely, 

N,ik + IT, g(0i), g(02)) = (1 - x{k'-'i<Pi - 02))) 

■{k + iT)^ai{k + iT, g(0i), g(02))e^('=+^^)l«(<^i)-«('^2)| 
with a,{k + IT, q{(t>i), g(02))) e 50(T2). 

Proof. We begin by analyzing the amplitudes of the Hankel functions. 
Lemma 4.4. There exist amplitudes ao, ai such that: 

(i) H^^\{k + it)z) = e*(^+^")^ao((A; + it)z), where (1 - x{k^~^ z))ao{{k + ir)z) G ^^^/^(M); 
{k+iT)H[^\{k+iT)z) = (A;+ir)^e*((^+*^)^ai((A;+ir)z), wzth {l-x{k^'^ z))ai{{k + iT)z) G 
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Proof. By the explicit formula (jHS)), 

We claim that the integral defines a polyhomogeneous symbol of order in (A; + ir)z as 
\{k + iT)z\ — »• oo. Indeed, applying the binomial theorem to the factor (1 — 2i{k+iT)z ^~^^'^ 
gives 
(49) 

= Ef=o ^Aik + tr)z)-^-^/' + R^{{k + it)z), 

where C(— l/2,j) are binomial coefficients, where cj are the resulting constants, and where 
RNis,{k + iT)z) = e~^s~^/'^RN{{k + iT)z,s)ds with -Riv((A; + iT)z,s) the A^th order 
remainder in the Taylor series expansion of (1 — 2i{k+iT)z ^~^^'^- evident that Rn{s, {k + 
iT)z) = 0{\{k + iT)z\~^^^^^'^^). Differentiating with k~^Dz similarly gives 



ik-^D)yoiik + tT)z) = J2j=o Cjiik + «r);2)-J-V2-l«l + {k-^DY,RN{k 



z 



(50) 

= 0(|fc2|-l/2+l«). 

Now 1 — x{k^^^ z) clearly belongs to ^^(IR). Since \kz\ > k^ on supp 1 — x{k^~^ z)i we 
conclude that 

{k-^DY^{l - x{k^'^ z))a{{k + ir)z) = 0(A;-(^/2+l°l^)), 

proving (i). 

By definition, H^\z) = —-j^Hq'^\z), so (ii) follows immediately from (i). □ 
We now complete the proof of the Proposition. The amplitude of A^^i is then 

(51) ai{k + ir, g(0i), g(02)) := ai{{k + ir)|g(0i) - g(02)|) cos^9i,2, 

where ■(9i,2 = /g(02) - g(0i), J^q{^2))- 

Since g : T — > dQ is smooth, the metrics |0i — 02|t and \q{4>i) — ^7(02)^2 are equivalent, 
and we use the former. Thus, we need to check that 

(52) \{k-'D^ra,{{k + zr)|g(0O - g(02)|) cos^i,2| < Ck-^'. 
By repeatedly differentiating (|THjl away from the diagonal, we find 

(53) l^?(g(0i) - g(02))| < c^lh - 02r-i"i 

It follows from (jlSp and by the chain rule that 
(54) 

{k-'D^ra,{{k + tTM<j>,)-q{<P2)\) = E,:,-<H(^"'^.))'«i((^ + ^^)^))|.=k(*0-.(*2)| 

X E,,...,,„|,|=|a|-.^«.,n^£=i(^~'^^)"^sin^i,2. 

By Lemma f4.4n i). 

(55) 1(1 - x{k'-\<Pi - <P2)m-'Dz)ya,i{k + 2r)z))|,=|,(<^,)_,(<^,)| | < C,k-^'. 
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Since \4>i — (p2\ > k on supp (1 — x{k^ ^{(f^i — 02))) it follows from that 

(56) 1(1 - x{k'~\(Pi - (P2)))ik-'D^r sin^i2| < q,|A;-l^^l^. 

The same kind of estimate is correct for the factor cos'i9i2. The proof folows from the 
combination of with the fact that (1 - xik^~ {(pi - ^2))) e ^"(T^). 

□ 

4.3. Layer potentials as semi-classical Fourier integral operators. In this section, we 
give a description of the layer potentials Si, T>i as Fourier integral operators from L'^{dQ) —>■ 
L'^{Q) which parallels that of Proposition 14.31 These operators also appear in the trace 
formula, and we will be needing the results of this section in §7. First, we must introduce 
suitable cutoff operators away from the diagonal. 

To define the cutoff operator, we need to introduce coordinates. We separate fl into two 
zones depending on the distance r{x, dQ) to dQ. We set 

( (dn), = {xen: r{x,dn) < e} 

(57) <^ 

[ a = n\{dn),. 

Here, e is sufficiently small so that the exponential map 

exp : N^{dQ) Q, expg(^^) rUg^^^^ = g(0o) + rUg(^cf>o) 

from the inward normal bundle along the boundary to the interior is a diffeomorphism from 
vectors of length < e to (dfl)^. Also, as above, i^q((f,o) is the interior unit normal at q{(po). 
The exponential map induces Fermi normal coordinates {(f), r) on the annulus (dQ)^, with 
r = r{x,dQ) the distance from x to dfl. We denote the Jacobian of the exponential map 
along the boundary by J{(f),r). We introduce the corresponding cutoff 

(58) xanix) ■= x{k'~'r), 

where it is understood that the cutoff is supported in a tube around the boundary where 
the Fermi normal coordinates are defined. We further introduce an angular cutoff as for the 
boundary integral operators. Near the boundary we may use Fermi normal coordinates (r, 6) 
and we introduce x{k^~^{(^ ~ 0) where q = q{(f)) in arclength coordinates. 

We then break up each layer potential into several pieces using the cutoffs. For instance, 
we first make a radial cutoff of the single layer: 

(59) x'9T\^))Si+ {I - x'aT'i^m^- 

We then further break up Xan'^^ i^))<Si into 
(60) 

xfn^\x))Si{k + IT, (r, e), g(0)) = xfn^\x))x{k'-' {9 - 0)5£(A; + ir, (r, 6), g(0)) 

+Xfn^\xm - x{k'-'{e - <P))Si{k + IT, (r, 9), q{<f>)) 

We do likewise with the double layer potential. 

We now show that, when suitably cutoff away from the diagonal singularities, the operators 
Si,Vi are semiclassical Fourier integral operators with phase function |x — g| : Qx dQ R. 
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The phase parametrizes the canonical relation 

(61) T = {(x,^^,q,-^^-T,)}cT*nxT*dn 

\x — q\ \x — q\ 

which is the graph of the interior-to-boundary billiard map (3, which takes an interior vector 
(x, |fE^) to the tangential component of the tangent vector(s) to the billiard ray it generates 

at its intersection point (s) with the boundary. The natural projections p : F — > T*Q, resp. 
g : r — > T*dQ have singularities at points where the billiard ray intersects the boundary 
tangentially. Away from these singular points, p resp. q are submersions with discrete, resp. 
1 dimensional fibers. Phases of this type arose in the work of Carleson-Sjolin (though not in 
relation to boundary value problems) and are discussed in detail in Sogge [So] . Since we are 
cutting off the diagonal we have eliminated grazing rays. In the case of convex domains, the 
cutoff layer potentials are non-degenerate Fourier integral operators, while in the non-convex 
they have degeneracies at points {x, |f5||!5')0) ^ T. These too will be eliminated when we 
microlocalize to the orbit 7. 
We have: 

Proposition 4.5. The operators 

• - XdQ^\x)))Ve{k + zt) (resp. Si{k + iT)); 

• X%T\x)){l-x{k^~\0-ct)))VE{k + iT,{r,e),q{ci))) (resp. Si) 

are semiclassical Fourier integral operator of order —1/4 ('Di), resp. —5/4 (Si), associated 
to the canonical relation P. More precisely, there exist amplitudes such that 

(t) {l-x'9n''\x)))Viik + iT,x,qm 

= (1 - Xan^\x))){k + zr)5Ai((A; + it), x, g(0))e^('=+*")l^-5WI, 

(n) X^n^\r))il-xik'-'i<p-9)))Viik + ^T,x,q{<P)) 

= {k + iT)^A2ik + ir, r, 9, 0)e^('=+*^)l'?(^)+''^'^w-«('^)l 

with Aj{k + ir, r, 6, (p) G ^^(T^). The analogous statements are true for Si. 

Proof. The proof is similar to that of Proposition 14.31 Due to the cutoffs, the kernel is in its 
semiclassical regime where the WKB approximation is valid and the orders of the amplitudes 
can be read off from Lemma lOl After substituting z = \x — q{(p)\ we get the amplitudes: 

• ao{{k + iT)\x — ^(0)1) for Si; 

• {k + iT)ai{{k + ir)|x — ^(0)1) cosZ(x — g(0), J-'q((i,)) for Vi. 

We note that the order convention on Fourier integral operators operating between spaces 
of unequal dimensions rii, n2 is that /c"-i/4+"2/4 |;ij2ies an amplitude of order on the stationary 
phase set defines an operator of order 0. Thus, the cutoff Vi has an order of 1/4 less than 
Ni. Since the normal derivative introduces a factor of {k + ir) in T>i but not in Si, the order 
of T>i is one higher than Si. 

To complete the proof, we need to check, that differentiations in k~^Dx, k~^D^ in case 
(i), resp. k~^D0, k~^Dr, k~^D^ in case (ii), lower the order of the amplitudes by k~^. The 
only difference to Proposition 14.31 is that ^q{4>2) — q{4>i) y ^q(4>2)) there is replaced by Z{x — 
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g(0), Uqifj,)) = '^{qid) + ruq(^ff) — q{(p), i^q{(i,))- As before, each derivative at most increases the 
number of factors of |x — g(0)|^^ by one. On the support of the cutoff, each such factor counts 
k^~^. Due to the accompanying factor of we see that each derivative k~^D decreases the 
order of the amphtude by k~^. 

□ 

4.4. Integral formulas. In §6, we are going to need an asymptotic formula and remainder 
estimate for some integrals involving Hankel functions. The reasons for the assumptions on 
the parameters a, b will be clarified at that time. The constraint l>5>l/2on/i = k~^~^^ 
originates in the following. 

Proposition 4.6. Let a g (-1, 1), let m>l,Qb>0 and let i? e N. Also, let x G C^(^) 
be an even cutoff' function, equal to one nearO. Then 

Xik-^x) cos{ax)H^'\bx)dx = + 0{k-^^\a - bf) + 0{k-^^), 

where the 0-symbol is uniform. The equation may be differentiated any number of times in 
a with the same remainder estimate. 

Proof. The result is suggested by the standard formula 



POO 

/ cos{ax)HQ^\bx)dx 
Jo 



Vb 



for the Fourier transform of a Hankel function [0 1 lAGj . To deal with the cutoff we proceed 
as follows. 

Since l[i oo](t)(i^ — l)^^^"^ is a tempered distribution on M, and since xik~^x) cos(ax) is a 
Schwartz function, we may replace the Hankel function by its Fourier integral formula (jH^ 
to obtain 
(62) 

Xik-^x) cos{ax)H^'^\bx)dx = J^{J^ xik'^x) cos{ax)e'^'''dx}{t^ - ly^'^dt 

= E± r r l[i,oo] - l)-'/'x{x)e^''^''^''^^dx. 

The phase has a (non-degenerate) critical point at x = 0, t = ±a/b (since 6 7^ 0). However, 
with our assumptions on (a, b), the critical point lies (slightly) outside of the support [1, 00] 
of the dt integral, hence the operator 

L = k^^-^—D^, 
bt ± a 

is well-defined on the support of the integral. We cannot integrate by parts in the dt integral 
due to singularities in the amplitude, but we can (and will) integrate by parts repeatedly in 
the dx integral with L, which reproduces the phase. 
The first time we integrate by parts with L, we obtain 

E± r r l[i,oo] - l)-'/\ix)^^D^e^''^'''^'^^^dxdt 

(63) = /„°°((6t ± a)-il[i,oo](t)(t^ - ly'/'dt 

+ E± r r i[i,oo](t)((t^ - i)-^/^(/^.x(^))^e^^^(^*±")^M. 
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We then integrate the second term by parts repeatedly with L and observe that no further 
boundary terms are picked up since [D^xi^)) = near x = 0. We then break up the dt 
integral into J^{---}dt + J^{- ■ - jdt. On the dt- integral over [2, oo], the factors of 
obtained from the R partial integrations render the integral absolutely convergent and the 
factors of show that it is of order 0{\k\~^^) with coefficients independent of (a, b). In the 
integral over [1, 2] we only have |6t±a| > |6 — a|. After R partial integrations we therefore 
obtain an estimate of k~^^\b — a\~^. for the integral. The same kind of argument applies if 
we first differentiate the integral n times in a, without changing the remainder estimate. 
To complete the proof, we note (with |AGj ) that 



(64) 



J2{bt±a)-\t^-ir/'dt=-^_ 



Vb^ - a? ■ 

□ 



Corollary 4.7. With the same assumptions as above, assume further that \a — b\ > 
Ck-^+^. Then 

J,^x{k-'x) cos{ax)H^^\bx)dx = + 0{k^('-''^), 

where the 0-symbol is uniform. Thus, for any 5 with 1 > 5 > 1/2, the remainder is 
rapidly decaying. The equation may he differentiated any number of times in a with the same 
remainder estimate. 

We will also need a slight extension of this result: 

Proposition 4.8. With the same notation and assumptions as above, we have 



Io° x{k ^x) cos{ax)HQ^\bVr'^ + u^) cos{au)dudx = — '^b^-a'^ + 0{k^^^ ^^), 
where the 0-symbol is uniform for a in the interval above. 

Proof. The proof proceeds much as before except that we now use the identity: (^, (14.16)- 
(14.18)) 

/ Hn^\bV + M^) cos( au)du = —i — , valid if b'^ — > 0. 

Details are left to the reader. 

□ 

4.5. The case of the unit disc. The only computable example is the unit disc, and the 
reader might wish to use it to check the calculations and estimaes. Due to the symmetry, 
N{k) is a convolution operator and there are simple commutation relations with the layer 
potentials. We have: 



(65) 



f Ar(A;,^,0) =kZnezHi'\k)j:^ik)e^-('"f^ 

Se{k, (r, e),<j)) = Hi'^ {k)Jn{kr) cos n{<f) - 6) 

[ Viik,ir,9),<P) =kY.^^^H^\k)J'^{kr)cosn{(l>~ 
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These layer potential identities follow from the so-called Graf addition theorem: 

(66) Go{k,x,y) = ^ H'^^\kmax{\x\, |?/|}) J„(A; min{|x|, \y\}) cos nZ{x,y), 

which is valid for general domains. 

5. MiCROLOCALIZING THE TRACE 

As discussed in the introduction, we are going to use (|7j) to calculate the wave invariants 
at Y- III Lemma f3. 31 we discussed the standard microlocalization of the trace to 7. We now 
prove that we can use a microlocal cutoff operator on the boundary as well as on the domain. 
This obvious sounding statement is not actually obvious in the present approach since the 
layer potentials are not standard semiclassical Fourier integral operators. 

Let us first review Lemma 13.31 in the language of layer potentials. It asserts that 

(67) p * Tr[V{k + iT){I + N{k))~^'''S{k + ir^ o (1 - x) ~ 0, 

where x{k) is a semiclassical cutoff in f2 to a microlocal neighborhood of 7 which has the 
form x(r, y, k~^Dy) near the boundary. Here, (r, y) are Fermi normal coordinates in Q'^ near 
the endponts of 7 , with r the distance to the boundary. Thus, differentiations are only in 
tangential directions. 

We now verify that the analogue of ()67|) remains correct if we use a microlocal cutoff on 
the boundary. The billiard map /? is a cross section of the billiard flow, and in this cross 
section a bouncing ball 7 corresponds to a periodic orbit which we denote ^7 G B*dQ of 
period 2. We choose the boundary parametrization so that = is one of the endpoints of 
the segment (which we also denote d'y) in Q. To microlocalize the boundary operators to the 
periodic point, we introduce a semiclassical pseudodifferential cutoff operator Xd^ 
with complete symbol Xd'y{(p,v) supported in := '■ |</>|, < e}. 

Proposition 5.1. Suppose as above that supp p n Lsp{Q) = {rL^}, and let xo-y 

be a semiclassical microlocal cutoff to the billiard-map orbit corresponding to 7. Then: 

p * Tr[V{k + iT){I + N{k))-^^''Sl{k + ir)**^ 

~ p * Tr[Ve{k + iT){I + N{k)y^'' o Xd^{k) o S{k + iTf. 
Proof. By ^ and by (inZ|)-(jSHI), we have: 

(t) Si{k + iTf = -fGoik + it) = /~ e''^''+'-^rn2Eo{t)dt, 

(u) lniVi{k + it) o (/ + Nik + it))-\x, q) = /J" e''^^+'-hnu2ED{t)dt, 

where r^j/M = dyu\da- We have subscripted the restriction operators to clarify which variables 
they operate on. The composition Vi o (/ + N{k))~^ o (1 — Xdj{k)) o Si^^ may be written 
(with the relevant value of k) as 

(68) / e^^C^+^^H / rn,2ED{t - s) o {I - Xd^){y, \Dt\-'Dy)rni o Eo{s)ds}dt. 
Jo Jo 
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We therefore have: 

p*Vio{I + N{k))-^ o (1 - Xd^{k)) o Sf"" 

(69) 

= ir rnu2EDit -s)o{I- xd^)iy, \Dt\-'Dy) o rf,iEo(s)rfs]e^('=+-)*rft. 

The statement of the Proposition is equivalent to: 

(70) WF[Tr [ rnu2ED{t-s)o{I~Xd^){yADt\-^Dy)ormEo{s)ds\n[rL^-e,rL^ + e]=^. 

Jo 

We now claim that the integrand 

V{s,t) = Tr rn,2ED{t - s) o (/ - Xd^){y, \Dt\-'Dy) o rm^o(s) 

is a smooth function for t G (rL^ — e, rL^ + e) and for s G (0, rL^ + e). Indeed, the singular 
support consists of (s, t) such that there exists a closed billiard orbit of length t outside a 
phase space neighborhood of Y which consists of a straight line segment of length s from 
a point X G r2 to a boundary point q, followed by a generalized billiard orbit of length 
t — s from q back to x. By our assumption on the length spectrum, the only possible orbit 
with length t G {rL^ — e,rL^ + e) is 7'' of length r; but the cutoff has removed this orbit. 
Since the integrand is smooth, the integral over s G [0,t] determines a smooth function of 
t & — e, rL^ + e). 

□ 

Remark 5.2. The Proposition is obvious in the case of the unit disc, although in this case 
it is only natural to cutoff in the frequency variable since all radial geodesies are bouncing 
ball orbits. Lemma VJ.'J^ and Proposition \5.1\ are then equivalent, since the cutoff to the right 
of St"^ has the form x{k~^D^) and commutes with Si^^ . 

6. Regularizing the boundary integrals 

The purpose of this section is to analyze the compositions in ()10|) with the semiclassical 
cutoff Xd-y on the boundary as semiclassical Fourier integral operators. Since the role of the 
imaginary part of the spectral parameter is not important here, we write it simply as r and 
only substitute r log k when it is needed in ^ 

For ease of notation, we write the terms of (fTUI) as 

(71) := N^ii) o Ar,(2) o . . . o N^t^M) 
and we put 

(72) |cr| = ^a^'^(0) = the number of A^^o factors occurring in No-. 

Proposition 6.1. (A) Suppose that is not of the form N^^ . Then for any integer 
R > 0, No o Xd-y{k + it) may be expressed as the sum 

No = Fo{k,^i,(P2) + KR, 

where Fo is a semiclassical Fourier integral kernel of order — |(t| associated to and a 

remainder Kr, which is a bounded kernel which is uniformly of order k~^. 

(B) Nq^ oxd-y ~ NoM oxd-y, where Nqm is a semiclassical pseudodifferential operator of order 
-M. 
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The proof will be broken up into a sequence of Lemmas. First we will consider the 
compositions NqoNi, NioNq without the cutoff Xd-y- Then we consider iterated compositions 
with the cutoff. Finally we discuss the special term Nq' o xa-y- 



6.1. The compositions NqoNi. Here we characterize the composition NqoNi. Essentially 
the same result holds for A^^i o A^g- 
The composed kernel equals 

No oN^{k + IT, 01, 02) := {k + it) xik-'^'{<Pi - 03))(1 - x(A:-^+'(02 - 03))) 

(73) 

Hi + ir)|g(03) - g(0i)|) cosZ(g(03) - g(0i), Ug(^^))Ni{{kfx + ir, g(02), g(03)). 

The somewhat technical nature of the following lemma is due in part to the lack of a cutoff 
to 7. 

Lemma 6.2. For any i? G N, there exists an amplitude AQi{k + ir, 0i, 02) such that 



• (i)NooNi{k+iT,(Pi,(l)2)) = (l-x(A;^"^(0i-02))A;(t-3'^)e*'=l^(<^i)-^('^2)Uoi(fc+^r,0i,02) + 
Kii{k + iT,(f)i,(f)2); 

• (11) Aoi{k + IT, (f),, 02) G 50(T2 X M). 

• (Hi) Kr is a bounded kernel which is uniformly of order k^^ . 

Proof. Using Lemma 14.41 and Proposition 14.31 we rewrite Ni in terms of phases and am- 
plitudes. The proof is then based on a change variables and on use of the explicit cosine 
transform of the Hankel function given in Proposition ()4.6|) to evaluate integrals involving 

the 'difficult factor' N^ik + zr), i.e. H^i\{k^i + zr)|g(0i) - g(03)|). 
It is convenient to first change variables 

03 ^ "i? = 01 - 03, 

and then to change variables with: 

|g(03)-g(0i)l, 0i>03 r |g(0i-^)-g(0i)l, ^>o 
-k(03)-g(0i)|, 0i<03 1 -|g(0i-^?)-g(0i)l, ^<o 



(74) u :-- 



In other words, we change from the intrinsic distance along dVl to chordal distance. Due 
to the cutoff, the variable u may be taken to range over [—k~^^^, k~^^^), so the change of 
variables is well-defined and smooth on the support of the integrand. The purpose of this 
change of variables is to simplify the difficult factor in (fTSj) : 

H^^\{k^^ + zr)|g(03) - g(0i)|) ^ {{k^^ + it)\u\). 
Now we consider the other factors. After the change of variables, we have: 
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(i) cos{Zq{(f)2) - q{(f)2 + ■&), ^q[<t>2)) \u\K{(f)2, u), with K smooth in u; 

{iii)Ni{{kfi + iT, g(02), g(0i - ^)) ^ e'^\i^'f'^'^-''^^^'^k''''"'A{k + ir, 0i, 02, u), 

where A/^ is a symbol in k of order 1/2 and smooth in u. 
Here, 

(75) a = sin^9i,2, with ^9i,2 = ^(g(02) - g(0i), t^qi<t>2))- 

These statements follow in a routine way from Proposition 14.31 and from the following 
identities (cf. jEP]): 

(z) g(0) - = (0 - - |/^(g(0))(0 - <l>')\i^') + - m 

{ii) |g(0) - = (0 - 0')' + O((0 - m 



(76) 



/N3 



(m) |g(0)-g(0')| = |0-0'|+O(0-0 



I (^^) ((g(0) - g(0')), ^qm) = - + O((0 - 0')=^). 



Here, T(0) denotes the unit tangent vector at g(0). It follows that K{u) = — l/2K(g(0j(,)) + 
(9(|Mp). For further details on the claims above, we refer to the Appendix |Z6j . 

Taking into account the factor (k + ir) in front of the integral, it follows that the composed 
kernel ()73|1 can be expressed in the form of Proposition ()(i.2|l (i), with 



1-5/ 



(77) 



Aik + IT, 02) = JZ Xik'-'u)il - xik 



X \u\e''"'''H[^\{k + iT)\u\)G{{k + ir), u, 0i, 02)rfu, 



01 -u)) 



where G{{k + ir), -u, 0i, 02) is a symbol in k of order 3/2 and smooth in u, and where 
x{k^~^u) = x{k^~^{(pi ~ 03))- The cutoff has been changed slightly under the change of 
variables, but for notational simplicity we retain the old notation for it. 

We now change variables again, u' = ku (and then drop the prime), to get 
(78) 



A{k+iT, 01, 02) = k~ 
with 6=1 + z(r/fc). 



x(rt)(l-x(A;i-^(02-0i-A;-iM)))|M|e*'^"i/fH6|M|)G(A;+zr,-,0i,</ 



Remark 6.3. We note that the entire amplitude now has order —1/2 at least formally. 
Since we have not yet microlocalized to 7, we will temporarily obtain a worse estimate, but 
in the final step we will show that this is the correct order. 

Since |m| < k^ on the support of the cutoff, we have ||| < k~^^^, and then the Taylor 
expansion of G{k + ir, u, 0i, 02) at m = produces an asymptotic series 

p 

fc-3/2 G'(fc + IT, p 01, 02) = J2 k-'^U^G^ik + IT, 01, 02) + k-PuPRp{k, p 0i, 02), 



n=0 
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where Gn{k + ir, 0i, ^2) is a symbol of order , and where Rp is the remainder, 

Rp{k, I, 01, 02) = -j\l- SY~^ G^P\k + IT, 01, 02)rfs. 

Since ^(A; + ir) is a symbol of order 0, G'^^\k + ir, s, 0i, 02) is uniformly bounded for |s| < 1 
and hence 

(79) |i?p(A;,p0i,02)| <Cp for|u|<fc^ 
Second, we Taylor expand the cutoff around m = to one order: 

(1 - x{k^-\^2 - 01 - k-\)) = (1 - x{k^-\<^2 - 0i)) + k-'uS,{k, p 01, 02), 

with 

u 

(80) ^i(fc,p0i,02) = ^' x{k^-\<Pi-<p2-s^))ds. 
We then write: 

Gil -x) = (Gp + R^m -x) + Sp) = Gpil -x) + GpSp + [i?p((l -x) + Sp)], 

where Gp is the pth Taylor polynomial of G. We claim that the first terms is a Fourier 
integral kernel of the type described in (i)-(iii); that the second kernel has this form but 
multiplied by cutoffs like x^(^^~'^(0i ~ 02)) which vanish except in a small band Ck-^+^ < 
101 - 02 1 < 10Cfc-i+^; and that the remaining two terms define an error of the type (iv). 
Consider the first term: 

poo 

J2 Gn{k + IT, 01, 02) (1 - x{k^-\<t>i - 02))) A;-"-2 / x{k-'u) \u\u^e'-^H? {h\u\)du. 

n=l -^-^ 

The integral may be expressed in the form 

(81) !Zx{k-'u)e'--H^^\h\u\)du\a=,^,,^^^^^^^ 

Since the terms in the finite part of the Taylor expansion of (1 — x(A;^^''(02 — 0i — k^^u)) 
vanish if |02 — 0i| < k~^^^ , we have a = sin'(9i_2 G (—1 + k^^^^, 1 — k^^^^). So we may apply 
Proposition ()4.6|) - Corollary ()4.7|) to evaluate (|HT|) asymptotically as 

ir X{kr^x) cos{ax)H^^\bx)dx\a=sin^^,2,b=a+iT/k) 



% 



^£L(i _ «2)-3/2 + o(A:-"^|(a - 1) + ^^|-(3+2")) 



= P„(a)(l - a2)-(3/2+n) + 0(A;-"^|(a - 1) + i^T^^^+^n)^ 

= (cos^9i,2)"(^+2") + 0(A;-"^A;(3+2n){i-5)). 

Here, P„ as an nth degree polynomial which we will not need to evaluate. 

Let us analyze the order in k of this part of the amplitude. We now put in the factor of 
^-n-2 front of ()81|) . We observe that cos'(?i,2 can be as small as k~^^^ on the support 
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of the cutoffs in the integral (which do not include a cutoff to 7). In view of the factor 
(cosi?i^2)~*-^^^"\ we estimate the contribution of the nth term to the amplitude as 

(82) A;~'^-2(cosi9i,2)-(=^+'") < A;-"-2A;(3+2n){i-5) ^ p(i-25)+(i-35)^ 

Each further derivative in k^^D gives a further factor of k~^. Since 6 > 1/2 the terms 
decrease in order with n, and since there are only a finite number of terms, we conclude that 
k~^/'^A{k + IT, 01, 02) lies in S'^^~^''\ This proves (i) - (ii) modulo the remainder estimate. 

Remark 6.4. In Lemma \6. 51 we will make use of the cutoff Xd'j o,t the end to eliminate 
these factors of cos 1^2- 

We next turn to the term GpSp, which is the most difficult of the remainder terms. As in 
(j8T|l of the previous step, the key point is to analyze the integrals 
(83) 

/o I-oo X{k~'u)x'{k'-'i<Pl - 02 - S^^)e^--H^'\b\u\)duds\a=sln,,^,,b=iWr/k). 

To deal with this integral, we need to adapt the method of proof of Proposition (|4.b|) - 
Corollary ()4.7|1 to take into account the additional cutoff. 
As before, we have: 

lo /o°° xif^~^ x)x' {(l^i — 02 — s|) cos{ax)HQ^\bx)dxds 

(84) 

= E± /o" /o /o" l[i,oo](t)(t' - l)-'/\ix)x'ik'-'i<t>i - 02) - sx)e^''(''^'^>dxdsdt. 

The phase is the same as in Proposition (j4.6|l . and we again integrate by parts in dx with 
L = k ^-g^D^. 

The first time we integrate by parts with L, we obtain 
(85) 

x'(A;i-^(0i - 02)) ± a)-H[i,oo](t)(t^ - l)-"'dt 

+ E± /o r r l[i.oo] - l)-''\D.M^W{k^-'{'i>i - 02) - sx)])^^e^^'i'^^->dxdtds. 

We then continue to integrate the second term by parts repeatedly with L. But unlike the 
case of Proposition ()4.6|) . we do pick up boundary terms each time when falls on the 
second cutoff factor. For instance, the next iterate produces the boundary term: 

fc- V'(A:^-'(0i - 02)) r((&t ± a)-2l[i,^](t)(t2 - ly^Ht 
(86) +k-' E± /o r r l[i,oo] - l)-"\Dl[x{x)^{k^-'{<P, - 02) - sx)]) 

{^fe''''^^'^''^^dxdtds. 

The boundary terms are similar to the kind in the GpX terms, except multiplied by x'{k^~^ 
and higher derivatives. This establishes the claimed form for the first two terms. 

As for the integral, we break up the dt integral into ^^{- ■ ■ }dt + {■ ■ ■ }dt as we did 
in Proposition ()4.6|) . The same observations show that after R partial integrations, the 
remainder term is bounded by k~^^\h — a\~^. for the integral. The boundary terms have the 
form of ^fj2_a2 times the series '^f=i k~^^ x^^Kk^~^ i'Pi ~ 02))- 
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We now sketch the estimate of the last two remainder terms. There are two kinds of terms, 
those from G which descend in steps of and those of x which descend in steps of k~^. If 
we retain the principal term of x and the first n terms of G, then we obtain a remainder of 

/oo 
x(rt)|u|u"i?„(u/A;,0i,02)e^""i/J^^(|n|)(iM 
■oo 

of the nth order remainder in the Taylor expansion of the amplitude. By fl79|l . we can bound 
the Rn and exponential factors by a uniform constant C. We can also bound the Hankel 
factor by j-ul^^/^ for \u\ > 1; the singularity at m = is cancelled by the m". Thus, the integral 
is bounded by k""-'"^ x"'~^^^'^dx ~ ]^-n-2^s{n+?,/2) ^ ^n(-i+5)-i/2^ ^ sufficiently large 

we have an arbitrarily small remainder. We can choose n to obtain the order stated in (iv). 

□ 

6.2. Proof of (A). By iterating the Lemma and composing with Xd-y we can improve Lemma 
l6.5l to the statement in (A) of the Proposition. First, some more notation. Any term A^. of 
(fTUI) can be expressed as a product 

(88) N'q^ o ■ ■ ■ N'q^ N{\ 

of blocks of Ao and Ai, with 'Yl!j=i — W\ ^'^'^ X]j=i('^i ~^ ^j) ~ The number r = r{a) 
counts the number of blocks of Aq. We now compose with Xd-y and successively eliminate the 
blocks Ag^ from right to left using one factor of Ai to the right. There is a slight notational 
problem since it could happen that ti = 0, in which case one should instead use the factor 
of Ai to the left in eliminating a block of Aq . Since the process is analogous in that case, 
we will assume for simplicity of notation that ti 7^ 0. 

Lemma 6.5. For any term A^- of mu\) except for Nq^ , and for any i? G N, there exist 
(si, . . . ,Sr) as above and amplitudes Aj{k + ir, 0i, 02), (with j = 1, . . . ,r) of order —Sj such 
that: 

• (i) N„ o xd^ = Nlr'^ oMr-o n'{-^~^ o ■ ■ . iV*^-i o Ml o N{^-^ o xa^ + K^r, where each 
Mj is a semiclassical Fourier integral kernel of the form 

Mj{k + IT, 01, 02) = k^'\l - x{k'^\4>2 - h))e"'^'^'^'^~''^^'^^A,{k + IT, 01, 02), 

with Aj{k + iT, 01, 02) e SJ'\T^ X M). 

• (iv) Kfi is a bounded kernel which is uniformly of order k~^. 

Proof. We work from right to left using the argument of Lemma 16.21 repeatedly to remove 
all of the Aq factors in each block. This can be done because we only used knowledge of the 
phase and of the order of the amplitude to obtain (i) - (iii) of the Lemma . In each such Mj 
we may choose R large enough so that the remainder for this block, when composed with 
the remaining factors of Aq, Ai satisfies (iii). We then define 

(89) F, = N'/-^ oMrO nI'-''^ o . . . iV{2-i o Ml o A^^^-^ o xa^ 
It is a semiclassical Fourier integral kernel with phase 

(90) /:<,(0i,...,0A/-|a|) = |g(0l) -g(02)| + ■■• + |g(0M-|a|) - g(0A/-k|-l) | • 

Now let us use the composition with Xd-y- Since A*'-~^oM^oA^*'-'"^o- ■ ■ A^""^ oMi o A^^^"^ 
is a semiclassical Fourier integral operator, its composition with Xd^y microlocalizes the kernel 
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to the periodic orbit of (3 corresponding to 7. That is, on the support of the cutoff Xd-^-, critical 
points correspond to Snell paths in which each hnk points roughly in the direction of 7. In 
calculating the order of the amplitude in Lemma [6. 2[ we had to take into account the factors 
of cos't?i_2, which we observed could be as small as k'^'^^ . However, upon microlocalizing to 
the periodic orbit, there is a uniform lower bound for cos'(9i^2; and hence these factors are 
bounded above. As remarked in the proof of Lemma I6.2t each removal of Nq decreases the 
order by 1 on the set where cos 191^2 is bounded uniformly below. Thus, the order of F^- is 
— |(t|. This improves the order estimate to the statement in (A) and completes the proof of 
this part of the Proposition. 

□ 

6.3. The term Nq^ o Xd'y- Our first step in proving (B) of the Proposition is: 
Lemma 6.6. A^^o o Xd'y is a semidassical pseudodifferential operator of order —1. 
Proof. The composed kernel equals 

A^o o xa,{k + IT, 01, 02) := {k + zr) x(fc-'+'(</>i " - x{k-'^\4>2 - h))) 

(91) 

Hi {{kfi + iT)\q{(l)3) - g(0i)|)cosZ(g(03) - g(</)i), 1/5(^3)) Xa7(^' ^('^a), g(03))#3, 

where 

X0,ik,q{(j)2),q{h)) = k f e'''^^^-^'^-^xa,{(f>2,v)dv 

Jr 

where Xd'y{(p,v) denotes (with a slight abuse of notation) the symbol of xo-y, a smooth 
phase space cutoff to (0, 0), /5(0, 0), the coordinates of the periodic point of period 2. 
We make the same change of variables as in Lemma f6. 21 which again takes 

' mi'\{kf, + zr)|g(03) - g(0i)|) - H['\{kfi + ZT)\u\y, 
< (a) cos(Zg(02) — ^(02 + ^q{<f>2)) ^ \u\K{(p2, u) , with K smooth in u; 

{iii)e^''^^^-^^)■'^Xa^{(t)2, v) ^ e*'=('^2-<^i)eifc^-r? ^ Xd^{<l)2, r])e'''^'^^-'t'^h'''''-'^A{k, u, ri)xd'y{<p2, r]), 

where A{k, u, rj) is an amplitude of order 0. 
Then fl91l ) becomes 

{k + tr){l - x{k-'+\<p2 - 0i)) 4 e^'('^2-'^i)a(fc + it, 0i, 02, ri)x^{<^2. v)dv, 

(92) 

with a(A; + ir, 01, 02, 77) = jr^x{k ^^^u)H'^\{kiJi + iT)\u\)\u\ e'^'^''^ Ai{k,u,7])dudr], 

for another amplitude of order 0. This is the same kind of integral we analyzed in (j87j) and 
(|8ip of Lemma [6.21 and we obtain the same description of a{k, 0i, 02) as for A{k + iT, 0i, 02) 
except that the parameter a in (j8ip is now rj. Thus, we obtain powers of (1 — if')~^, which 
blow up on the unit sphere bundle of dVL. These points of course correspond to tangential (or 
grazing) rays, and on the support of the cutoff — is uniformly bounded. Thus, we 

obtain the statement of the Lemma precisely as in Lemma f6. 21 and its improvement Lemma 
I6.5l with the cutoff. The change of variables eliminated two powers of k, leaving an amplitude 
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of order 0; an amplitude of order defines a semiclassical pseudodifferential operator of order 
-1. 

□ 

To complete the proof of B it suffices to iterate Lemma Ffj. 61 On each application, we have 
a new cutoff operator, but there is no essential change in the argument. This completes the 
proof of Proposition 16.21 

7. Regularizing the interior integral 

We further need regularize the integrals involving the outer factors of Vi{k + it) and of 
Si{k + irY^. Since we are taking a trace, we can (and will) cycle the factor of Vi{k + ir) to 
the right of Sf^{k + iT)x{k) to obtain an operator 
(93) 

Sex{k)Ve: L'^{dn) ^ L^{dn), Stx{k)Vl{q,q') = I Go{X,q,x)xik-^D.^,x)d^Go{X,x,q')dx. 

Jn 

on dQ (or S^, after parametrizing it). Here, x{k) is a semiclassical cutoff to a neighborhood 
of a periodic orbit 7. For simplicity we will assume that 7 is a bouncing ball orbit, although 
the same method would apply to a general periodic reflecting ray. 

We note that use of the inside/outside duality as in |Z5j would in effect make this section 
unnecessary, at the expense of forcing a stronger hypothesis on the simplicity of the length 
spectrum. Indeed, the integrals over the inside/outside would add up to the integral in 
but over instead of fl. This integral is easily evaluated to be N' [k + ir). We will now show 
that the integral over (with the cutoff in place) produces a similar kind of semi-classical 
Fourier integral operator. 

Before stating the precise results, we give some heuristics on the composition (j93p . In 
Proposition 14. 5( we described the layer kernels away from their diagonal singularities. We 
are now including the latter singularities. In addition to F of (|UT|. the wave front description 
now also includes the relation 

(94) A, = {(g, g, v) ■■ i\dn = r]} C T*n x T*9fi, 

which carries the singularities of the kernel. Intuitively, iS£, T>^ are singular Fourier integral 
kernels associated to the union (which we write as a sum) of the two canonical relations 
As + F, and hence the composition St o x{k) ° T^d- should be associated to the composition 

(A, + F)* o (A, + F) = A* o A, + A* o F + F* o A, + F* o F. 

It is clear that 

F* o F = A* o F = F* o A, = F^, 
the graph of the billiard map. Also, 

A* o A, = Ae, 

the diagonal oiT*dVt. We therefore expect the composition to contain these two components, 
precisely as N'{k + ir) does in the combined inside/outside case. The following proposition 
confirms this. It also produces a cutoff function, which confirms Proposition 15. II 
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Proposition 7.1. Six{k)'Di ~ Xd'y{kD0,9)[Do + Di], where Di is a semiclassical Fourier 
integral operator of order —1 associated to the billiard map, where Dq is a semiclassical 
pseudodifferential operator of order —3, and where Xd-y{kD0,6) is a microlocal cutoff to the 
(3-orhit o/(0,0). 

Proof. We break up each term as in ()59|) and then ()60|) and analyze each one separately. 
7.0.1. The most regular terms. The regular terms are those of the form 

(95) / {I - Xm^\x))fG,{X,q,x)x{k-^D^,x)d,Go{X,x,q')dx. 
Jn 

or where the cutoff has the form (1 — Xm^\x)))Xdn^\x))- 

Lemma 7.2. The regular terms of the form / l.9,5j) define Fourier integral kernels of the form 

where A is a semiclassical amplitude of order —1, and where Xd-y is a cutoff to 7. Thus, / fP3]) 
defines a semiclassical Fourier integral operator of order —1 associated to (3. 

Proof. Each factor of the product is described by Proposition 14.51 Further, we may take to 
be the product of a frequency cutoff x{k~^Dg) and a spatial cutoff ip-yix) to a strip around 
7. Thus, the integral has the form: 

(96) {k + zr) / (1 - ;^|^^+\a;)))2e^(('=+-)[l-<?l+l'''-IV,(x)x(^^)A(fc, x, q, q')dx 

where A is a semiclassical smooth amplitude of order —1. The stationary phase set is defined 
by 

/ ^ I I 7 1 'I ^ ~ ^ 2; - g' 

\x : dx\x — q\ = —dx\x — q \ <^=^ 1 = —- r|. 

\x — q\ \x ~ q'\ 

First, we see that x G qq' (the line segment between q, q'). The ray qq' is constrained by the 
cutoff to point in the nearly vertical direction in the strip containing 7, hence \q — q'\>C>0 
on the support of the cutoff. (Note that the phase is the sum, not the difference, of the 
distances since we are taking the transpose, not the adjoint, of Si{k + it). For the adjoint, 
the critical point equation would force q = q'). 

The ray qq' is thus a critical manifold of the phase, and the phase equals |g — g'| along it. 
To prove Lemma 17.21 we show that this critical manifold is non-degenerate and determine 
the amplitude by stationary phase. We choose rectangular coordinates (s, t) oriented so that 
the t axis is the ray qq' and so that the s-axis is orthogonal to it. Then q = (0, b),q' = (0, b') 
for some b,b' eR and the phase may be written \E' = (s^ + (t - 6)^)^/^ + (s^ + {t - b'^y/'^. 
Hence, = s[(s^ + {t - b)"^)'^^"^ + (s^ + (t - 6')^)""^^^] and on the stationary phase set s = 
the second derivative is simply \I^"^(0,t) = [((t - bY)~'^/^ + {{t - 6')^)"^/^]. It is clear that 
^^"^(0, t) is uniformly bounded below and since |6 — 6'| > C on the support of the cutoff it is 
also uniformly bounded above. Stationary phase introduces a factor of A;"^/^ so the resulting 
amplitude is a product of k^^'^ with an amplitude of order —1 and hence the composition 
defines a semiclassical Fourier integral operator of order —1. 

□ 
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7.0.2. Singular terms. Therefore we only need to consider the integral 
(97) / [xfn^\x))]^Go{X, q, x)x{k-'D,, x)9.Go(A, x, q')dx. 

JQ 

With no essential loss of generality, we redefined the cutoff to remove the square. We write 
X = (r, 6) and then break up the integral into the sum of four terms corresponding to the 
cutoffs: 

• (i) xrn"\r))x{k'-'{9 - ct>)x{k'-\e - ^ 

• (ii) X'on^'irmi - x{k'-'{9 - <P))x{k'-\e - <P')- 

. (iii) xfn^\r))x{k'-'{e - m - xik'-\e - 0')); 

• (iv) xrn"\r)){l - x{k'-'{e - 0))(1 - x{k'''{e - 0'); 

7.0.3. Codimension zero case. This refers to case (iv). There are no singularities in the 
integrand due to the cutoff. We can use the WKB approximation in each Green's function 
and apply the cutoff operator in the smooth variables 0, (f)' to obtain integrals of the form: 

/; C X%T\rm - x{k'-\e - mi - X{k'-'{e - 0')) 

(98) 

x(V^|x(r, e) - g')ei({fc+i^)[I^M)-g|+k'-^{r,0)|]^(^^ q')rdrde. 

There is no essential difference to the regular terms in (|95|). We therefore have the same 
result. 

Lemma 7.3. Integral i98^} defines a Fourier integral kernel of order —I of the same form as 
Lemma \7.^ . 

7.0.4. Codimension one case. This applies to cases (ii) - (iii), which are quite similar although 
not identical. We do case (ii); case (iii) is similar. 
We have: 

lo C xfu^'imi - x{k'-'{e - <p'))x{k'-'{e - 0) 

(99) (n) = 

Go(A, g(0), r, e)x{k-'D^, x)d,Go{\ (r, 6), q{(t)'))rdrde 
Lemma 7.4. Integrals \99\} define Fourier integral kernels of the type 

k''\a,{q, g', ^^)e^('=+-)I^W-^(<^')lA(fc, r, g(0), g(0')), 

with A an amplitude of order —2, i.e. they are semiclassical Fourier integral operators of 
order —2 associated to (3. 

Proof. The proof is reminiscent of that of Lemma but is somewhat more complicated. For 
the sake of brevity, we concentrate on the new details and do not discuss the error estimate, 
which is similar to that in the proof of Lemma 16. 2|) . 

We substitute the WKB approximation for dyGo^X, {r,9),q{(f)')) (see Proposition 14. 5p but 
not for Go{X,q{(f)),r,6) and apply the cutoff operator to the WKB expression. The integral 
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of concern is thus: 

(Unni) := {k + tr) /; xik'-'r)x{k'-'iq{<Po) + rUgiM ' 9(0i))) 

(100) 4'\ik + ^r)|g(0) + rz.,(,) - q{9)\) xi^§^SS^\ " ^'^^ 

A{k + IT, e, r, 0')e^('=+*")l''W+'-^.w-''('^'))l(l - x{k^'^\q{0) - q{(p')\)rdedr, 

where A G SJ^^^ 

We cannot use the WKB expression for the factor H^^\{k + iT)n\q{6) + rz/q(0) — q{(t))\), so 
we deal with it by changing variables and explicitly integrating. The change of variables is 
given by (r, 0) — >• (r, m), with : 



(101) u{e) :-- 



(|g(^) + rz/,(,)-g(0)p-r2)V2, 
-{\q{0)+rv,i^e)-q{<l>)?-r^f'\ & > <l> 



Here, (p — 6 and -u range only over {—k^^^^,k^^^^). We claim that m(0) is smooth and 
invertible with uniform bounds on derivatives as r varies on [0, eo]. 

For the sake of brevity, we only verify this in the basic case of a circle of radius a and refer 
to |Z6j for the routine extension to the general case, which only requires showing that the 
quadratic approximation is sufficient to determine the smoothness of the change of variables. 
We recall here that we have microlocalized the integral to 7 and that dfl has non-vanishing 
curvature at the reflection points of 7. We then have: 

{\qi9) + ru,^e) - q{<P)? - r') = 2(1 - r)(l - cos(0 - 9)). 

For 1^ — 01 < —k~^'^^, (1 — cos(^ — 0)) has a smooth square root, given by the stated formula. 

We now consider the effect of this change of variables on the remaining factors. For the 
exponential factor, we have 

I gi(fc+jr)|g(e)+ri/,(e)-g(9i'))l ^ik\q{(f>)-g{(p')\^ik{ua+rVl-a^) 0^ 

where a = sin Z(g(0) — g(0'), i^q((f,)) and where Ai{k, /i, 0, 0', u, r) is a polyhomogeneous sym- 
bol of order in k. We note that F(0, 0) = 1,G(0,0) = -l/2fi;(g(0i)). Also, the cutoff 
transforms as: 

X87 I (a\ I TJiV\ ' ^ ^9'<^T7T\ TJiV\ ' ^s))K{k, r, q,(l),(f) ,6), 

where K{k, r, q, (p, 0', 6) is a smooth semiclassical amplitude whose leading order term equals 
1. 

After multiplying together these amplitudes, and rescaling the variables u ku,r kr. 
In the regime (p — 6 = 0{k~^~^^), we obtain we obtain an amplitude A{k,(j),(j)') such that: 
(CnnD = Xd^iq, q\ ^)e^'=l'?W-''(^')lA(fc, 0, 0'); with 

A{k, 0, 0') = fc-i x('t"^«)x('^"^^)e^^^""+"^^^ 

(102) ^ 

H''Q\bVu^^)Ai{{k + ir), M, 0, 0', r, fi)))dudr, 
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where b = + ir/k) and where Ai is an amphtude of order —1/2. As a check on the 
order, we note that the normal derivative contributed a factor of {k + ir) and the change of 
variables put in leading to the stated power of k. 

— 1/2 

To complete the proof, we need to show that A & Sg (T x R). We analyse this integral 
()102|) . working by induction on the Taylor expansions of Ai in the integral 



OO POO 



We recall from Proposition ()4.8|1 that: 



J{r; a, b) := / H'^\h\/r'^ + u^) cos{au)du = —i — , valid if b'^ — > 0. 

Jo \Jb^ - 

Taylor expanding the cos Zg(6') + rVf^[Q^ — Vq^^)) leads to the following integrals: 

As in Proposition ()4.6|) . the cutoff factor gives lower order terms in k. 
We then integrate in dr. The basic integrals are : 



(103) 



To leadmg order m k the exponentials e ^ and cancel, and we find the integrals 

grow at the rates A;~^+'' cos-i?^/^. By differentiating in (01,02), can obtain any term in the 
Taylor expansion of F with remainder estimate. Putting in the higher order terms in the 
Taylor expansions of F just adds lower order terms in producing a symbol expansion 
as in the previous cases. Due to the cutoff factor X07( |g(^)lg(0'j| ' ^e)); the factors of cos'(?^\ 
are bounded above. The details are now similar to those in the proof of Lemma (6.21 

□ 



7.0.5. Codimension two. The final integral we must consider is 



(104) 

/o /51 xin^\x))x{k-^+'{9 - mx{k-'^'{0 - 0O))G'o(A, q, x)x{k-^D,, x)d,G,{\ x, q')rdrd9. 

In this case, we cannot use the WKB expansion for either Green's function and must use 
integral formulas for products of Hankel functions. This is the most complicated case, and 
it is the one producing the pseudodifferential term in Proposition 17.11 
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Lemma 7.5. Integral ^1U4\ ) defines a semiclassical pseudodifferential operator of order —3 
with kernel of the form a{(j), D^) o Xd-y, ^-e. with kernel of the form 

where A is a semiclassical amplitude of order —3. 

Proof. The integral we are considering is: 
(105) 

(Uni = (fc + tr) /; x(A:^-V)x(A;i"'5(g(0) + ru^^e) - g(0))) 

Xik'"'{qie) + ru^^e) - g(0'))) + ^r)\q{e) + ru^^e) - g(0)|) 

x{k~^De, e, r)H['\{k + zrMd) + rv.i^s) - I) ^9 (^) + rv^^e) - q{<f)'), v^^e))rdedr. 
We substitute the Fourier integral formula for x{k~^Dg,9,r) to obtain: 

(Uni := ik + tr) j; x{k'-'r)x{k^-\q{cp) + ru.^e) - g(0))) 

Xik^-'iqi9) + rz/,(,) - g(0'))) i^i'Hl^ + ^rMO) + ru.^e) - 

(106) 

x{p0,e,r)e'^P<^('-''m^^\{k + ir)\q{d')+rv,^e')-qm\) 

Zq{6') + rug(^g>) - q{(f)'), Ug(^e'))rd6drd6'dpe. 

We then make the change of variables {9, 9', r) — > {u, u', r) defined in ()107|) with respect 
to the pairs {9;(f),r) {u,r) and {9';(j)',r) — > {u',r). Thus, 

r (|g(^) + rz/,(,)-g(0)p-r2)V2, 

(107) M := < 

I -(|g(^)+rz/,(,)-g(0)p-r2)V2, 

while m' is defined similarly with [9', 0') replacing {9, (p). We note that m ~ (^ — 0) since 
= m - qm' - 2rK{<PM9) - qm9 - <P) . sin Z ^^1^1 ~ ^1^1^ , ^.,(,)) + • ■ • 
This gives 

m = j: jz iZo iz xik'-'r)x{k'-'u)x{k'-'u') 

X{pe, 9{u, r, 0)^ r)e^''P»(^("'<^''')-^'("'''^'''')) 
H[^\{k + iT)Vu^ + r^)Hl^^\{k + tT)^{u'y + r^) 

C^^''''T,:ZT''^^ )B,{u, u', r, 0, <P')dudu'drdpe. 

Here, 5^ is a smooth, polyhomogeneous amplitude. We recall that x{.P9^u{9,r, (l)),r) is 
compactly supported in pg so that there is no problem of convergence of the dpo integral. 



(108) 
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We then change to scaled variables k~^u, k~^u', k~^r and expand 

(109) 9{k-%,k-^r,^) = k-% + ^ + --- 
to obtain 

(110) (IinSD= x(P9,O,0))e^'=f«(^-<^')Afc(0,0',pe), 
where 

A,(0, <P',pe) = k-' /; JZ x{k-'p)x{k-'u)x{k-'u') 
X'ipe, 0{k~^u*, k^^r*, (p), fc-V)e^P«("-"') 

Bk{k^^u, k^^r, (p, (f)', k~^u')dudu'dr. 

We would like to prove that Ak is a symbol. As above, we work by induction on the Taylor 
expansion of F, G, B^. Polynomials in (u, m', r) may be expressed as sums of derivatives with 
respect to parameters (oq, oi, ho, hi) at = ai = of the integrals 

k-' lo IZ IZ X{k-'p)x{k-'u)x{k-'u') 
X{pe, 0{k-\ k-^r, 0), A;-V)e'P»("-"'^ 



or the same with (- 



replacing 



We now do the du\ and du^ integrals first as in the 0o ~ 0i case. The zeroth order terms 
in the usual Taylor expansions produce dr integrals of a form similar to (jlOHj) : 

' ^^^/~e*'^2r^(fc-V)e-^"[V^rfr 



The rest proceeds as in the proof of Lemma 16.21 We conclude that 



□ 



8. Tail estimate 

In this section we provide the remainder estimate. As discussed in the Introduction, the 
remainder comes from the Neumann series: 
(112) 

Mo 

(/ + iV(fc + 2r))-i = ^(-l)^iV(fc + zr)^ + 7^Afo, 7^Mo = Ar(fc + zr)^°+^ (/ + iV(fc+2r))-^ 

M=0 
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It is now important to take the imaginary part to be of the form r log k. We further specify 
p to have the following properties: p{t) = po{t — L) where po G C^(M) is non-negative and 
supported in an e-interval around 0. 

Lemma 8.1. For any R, there exists Mq = Mq{R) and 

Trp * Ve{k + IT log A;)7^Mo + log k)Ve{k + ir log kf^xik) = 0{k~^). 

Proof. By Proposition 15. ![ we can insert the cutoff Xd^y into the trace to obtain: 

Trlnp * Vi{k + iTlogk){I + N{k + iTlogk)-^N^^''Xdy{k) 

(113) 

o(fc + ir log k)Si^{k + ir log k)x{k) 

where S^{k,q,x) = r^Go{k,x,y) is the transpose of the single layer potential. The cutoffs 
may be inserted at various points, and for simplicity of notation we often write only one in 
(the other can then also be inserted). 

We regard the trace as a Hilbert-Schmidt inner product for Hilbert- Schmidt operators 
from L'^{dQ) to L'^{Q). For ease of notation, we the remainder at step Mq + 1 and the trace 
is then 
(114) 

p * {Vi{k+iT log k) {I +N{k+iT log k)-^N{k+iT log k), [N^° {k+ir log k)xd^Sf{k+iT log k)]*)Hs- 

Since |p| = po is a probability measure, we can estimate the convolution integral by Cauchy- 

Schwarz inequality as 

(115) 

Po * \{Vi{k+iTlogk){I+N{k+iTlogk)~^N{k+iTlogk), [N^'^°{k+iTlogk)xd-ySt{k+iTlogk)]*)Hs 
We further apply Schwartz' inequality to the inner product to obtain the upper bound 
Po * \\Vi{k + iTlogk){I + N{k + iTlogk)-^N{k + iTlogk)\\jjs 

(116) 

■ I liV^^o (A; + log k)xd,St{k + IT log k)x{k) Wis- 

We now separately estimate each factor. The first estimate is crude but sufficient for our 
purposes. 

Lemma 8.2. \\Vi{k + iTlogk){I + N{k + irlogky^NHns = O(^). 

Proof. We first use the inequality IIAfillj:/^ < | |y4| 1 1 where H^H is the operator norm 
to bound 

\\Vi{k+iTlogk){I+N{k+iTlogk)-^N\\Hs < \\Vi{k+iTlogk){I+N{k+iTlogk)-^\\ \\N\\hs- 

Here, || ■ || denotes the L'^{dQ) L^(fi) operator norm. 
By Proposition 1)4.11) (ii), we have we have 

\\N\k + tTlogk)\\Hs = Oik'/''). 

By fl37|) - ()HHj) the norm \\Vi{k + ir log k) {I + N{k + IT log k)-^\\ is the norm of the Poisson 
operator. We claim that 

(117) \\PI{k + iT)\\man)^mn) < r-^k"^+\ 

This estimate is probably crude but it is sufficient for our purposes. We are not aware of 
any prior estimates on mapping norms of the Poisson kernel in the k aspect. 
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We first rewrite the Poisson integral of ()H7|1 as 

(118) Plik, X, q) = rnX^Gnik, x, y) 

where Xi, is any smooth vector filed on VL which agrees with on dVL. The operator r^Xy 
is (roughly) of order 3/2, so we rewrite the composition as 

(119) PI{k + IT, X, q) = rnX^A-^^^-'A^^^'Gnik + it, x, y). 

Here, Aq is the Dirichlet Laplacian and the fractional powers are defined by the method of 
Seeley jSl]. Both Xy and operate in the y variable. 
To prove pi7j) . we use the (well known) fact that 

is a bounded operator. Furthermore, 

is bounded for any e > 0. (See |LMj . Theorem 9.4 for proof of the continuity for e > 
and Theorem 9.5 for proof of lack of continuity if e = 0). So the first factor is bounded 
independently of + it. 

For the second, we use that Ali^^'Gn{k + iT,x,y) is the kernel of A^''^^^(An + (A; + zr)^) ^. 
This is clearly a bounded normal operator on L'^{Q), so its L^(fi) —>■ L'^{Q) operator norm 
is given by 

(120) \\A'J'+\An + {k + tT)')-'\\L^n)^LHn) = "^f^ J^^^fjT^-^^- 

It is elementary to maximise this function, and one finds that the order of magntitude of the 
maximum occurs when Aj ~ /c, and it then has the form 

(121) 

|rfc| r 

This completes the proof of the Claim and hence of the Lemma. 



□ 



We now give the crucial estimate. It explains why we did not need sharp estimates in the 
previous step. 

Proposition 8.3. For any R, there exists Mq such that: 

\\p * N''^^{k + tT\ogk)x9,S£\k + tThgk)xik)\\ls = Oik-''). 

Proof. □ 

We estimate \\Si{k + ir log k)xik-'^D^)\\ by a power of k as above, and thus reduce the 
proposition to estimating 

(122) \\\p*N{k + tT\ogk)''''xdj\\HS 

where the HS norm is now on Hilbert-Schmidt operators on L'^{dQ). We write out the 
Hilbert- Schmidt norm square as the trace: 

(123) Trp * N{k + IT log k)^^''Xd^X*d^N{k + ir log k)*^° 
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where the trace is on L'^{dQ). 
By Proposition 16.11 

(124) Trp*F,,X8,xl,F:, 

(Tl,(72:ZMf,-+{0,l} 

plus errors which may be assumed to be 0{k^^). We recall that is a Fourier integral 
operator of order —\a\ associated to p^^o-W\^ jj^ addition, the phase of F^j has the form 

(125) ik{Ca-, - €^2) - rlogk{C^-, + C^^), 
where 

■ ■ ■ ,'?Afo-k|) = kl - Q'2| H \- I^Mo-lal-l - Q'A/o~|a|| 

is the length of an Mq — |cr|-link. The sign difference in the first term reflects the composition 
of A^*^" with its adjoint, and since the second term comes from the real part of the phase 
there is no sign change. 

We now estimate the traces by applying apply stationary phase. We are only interested 
in the order of the trace and not in the coefficients, so we argue qualitatively. The terms 
of the stationary phase expansion of Trp * -^0-1X97X37-^^2 correspond to the fixed points of 
the canonical transformation j3^-^o-Wi\ q ^-(A/o-|cr2|^ _ ^-|cri|+|o-2| -^^j^ich lie in the support 
of the cutoffs Xd-y- By assumption, the only periodic point of f3 in the support of the 
cutoff is the period 2 orbit corresponding to 7. Thus, other critical points (which we will 
call general critical points) can only occur when \ai\ = \a2\, which we henceforth write as 
\a. Equivalently, the general critical points correspond to a closed path obtained by first 
following any Mq — |cr|-link Snell path from a variable point qi G dfl to an endpoint q^, 
and then reversing along the same path back to qi. The general critical points thus form a 
non-degenerate critical manifold parametrized by {qi,qo) G dfl x dQ. 

Now consider critical points in the support of X97- It vanishes unless the first link points 
roughly in the direction of the first link of 7. Since critical paths are Snell, this forces all 
links to point roughly in the directions of links of 7. It follows that all links of critical paths 
(including Y ^ind the general ones) have lengths ~ CL^ for some absolute constant C > 0. 
Note that C = 1/2 for 7, so this is approximately correct for all links. Thus, the imaginary 
part of the phase introduces the damping factors 

(^126) g-2Crlogfc(Mo-|o-|)L^ 

into the stationary phase expansion. It follows that 

(127) Trp * F,,Xd,Xh^F:^ ~ ^-2C.(Afo-k|)L,^-2|H^-(Afo-kl)+i j « 

for some smooth density a. We arrived at this order due to : 

• the order — |cr| of F^- fProposition 16. 1|) : 

• The damping factor p26|) : 

• Application of stationary phase to an 2Mq — (|cri| + |cr2)-fold integral with a 2- 
dimensional non-degenerate stationary phase manifold. 
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Given i?, we need to choose (Mq, r) so that 

-2Cr(Mo - \a\)L^ - 2\a\ - (Mq - \a\) + 1 < -i? 

^ (Mo - \a\){2CTL^ + 1) + 2|(t| >/? + !. 

for every a. For the case \a\ = Mq, it suffices to have Mq > (i?+l)/2. Otherwise, Mq — |cr| > 1 
and it suffices to pick r > With these choices of (Mq, r) the inequahty is true for all a. 

□ 

□ 

9. Completion of the proof of Theorem 11.11 

To complete the proof it suffices to determine the trace asymptotics of the regularized 
finite sums 

Em=o E<x:{i,...,Af}^{o,i} Trp * Ar,(i) o Ar,(2) o . . . O N^^M) o Xd^ 

(128) 

oSi^''{k + iT log k)ox{k)oV£{k + iT log A;). 
We recall (see^) that p * A{k) is short for 

(129) / p{k - p){p + iT)A{p)dp. 

Jr 

Proposition 9.1. Let be a periodic m-link reflecting ray, and let Y be its rth iterate. Let 
p G C^irLry — e, rL^ + e) he a smooth cutoff, equal to one near rL.y and containing no other 
lengths in its support. Then there exist coefficients a^±r ^.j such that 

Tr p * iV^d) o iV,(2) o . . . o N„^M)Xa'ySf'-{k + it) o x{k) o ve{k + it) 

g(ifc-r log fc)rL, k-^'^^i^f^^ia^r^^j + a^-r^,j)k'^ + 0{k~^)}, (M > mr), . 

Onikr^), M <mr 

Here, \cr\ = {j : a{j) = 0}. In the special case where M = mr and a{j) = 1 for all j , we 
write a^±r o-j) = bj±r.j. 

The proof consists of a sequence of Lemmas. We begin by collecting the results of Propo- 
sitions 16.11 and 17.11 

Lemma 9.2. Tr p* N^^i^ o N^^2) o o ■ ■ ■ o N^^M)Xd',<St'{k + irhgk) o xik) o {k)V£{k + it) 
is a finite sum of oscillatory integrals of the form 

UM-\a\+?,)l2 r r r ifc[(l-/.)t+M£<.(g(</.i),...,q(</.M-i<.|)]p-^logfc£(g(0i),...,g(</.M-|<.|)) 



X{(l{(l)i) - ?(02), 0i)v4(A;/i, 01, ... , (t)M-\a\)p{t)dtdiid(t)i ■ ■ ■ d(j)M-\a\, 

where 

^a{q{(j)i), • • • , g(0M-|a|) = |g(0i) " qih) | + " " " + |g(0M-|a| " g(0i) I 



X(9(0l) -g(02),01 
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and where A{k, 0i, . . . , (j)M-\a\) ^ 'S'^''^'- 

Proof. There are two somewhat different cases, namely the case where |cr| < M and the case 
where = M. 

9.0.6. Case (i): \a\ < M. By Proposition 

(130) Ar,(i) o N^(2) o ■ ■ ■ o N„(^M){k, 01, 02) = F^{k, 01, 02), 
where is a semiclassical Fourier integral kernel of the form 

(131) F„{k, 01, 02) = e^('=+-)l9(*O-9(02)|^(^i-5(^^ _ 02))A^(A;, 0i, 02), 

with A„ E /^l°"l(T^). By Proposition 17. H the full composition 

(132) iV^d) o iV,(2) o . . . o N^^M)Sf''{k + tT)xdj o xik) o P£(fc + zr) 
has the form 

(133) F^{k)oxdjo[Do + Dil 

where the operators Dq,Di are from Proposition 17.11 Di puts in an extra factor of T and 
an amplitude of order —1. All of these operators are semiclassical Fourier integral operators 
with symbols in 5*1 and therefore they can be composed in the standard way. Therefore we 
have: 

(134) F^{k) o xa^Si''-{k + iT log k)ox{k)oVi{k + iT log k) G /^"'"'"^(T^). 

We then unravel p* and recall that the factor {k — n) raises the order by one. Finally we 
change variables fi ^ kp which again raises the order by one. 

9.0.7. Case (ii).- |cr| = M. By Proposition Nq'^ o xa-y is a — Mth order semiclassical 
pseudodifFerential operator, so it suffices to consider composition of the form A_mo Xd'y[Do + 
Di]. The statement is clear in this case. 

□ 

The next step is to show that the stationary phase method applies to oscillatory integrals 
in Il^{Y, $). This is almost obvious, but for the sake of completeness we include the proof. 

Lemma 9.3. Let /^(a, $) G Ig^{Y,^), and let C$ denote the set of critical points of ^, 
and assume that C$ is a non-degenerate critical manifold of codimension q. Suppose that 
1/2 < 5 < 1. Then: 

/fc(a,<l>) ~ I 

[ RemM,Ria, $, R)k-^, C$ = 

where aj is a polynomial in the jet of the amplitude a and phase $ at C$, and where the 
remainder Rem ji{a,^, k) < | |a| |c3jj + 1 1$| |(73ii. 

Proof. We consider an oscillatory integral 
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with $ G S^, with A{k,y) G S^{Y). By assumption, the critical set = {y : Vy$ = 0} 
of the phase is a non-degenerate critical manifold. We choose a cutoff ijj supported near 
and write 

I{a, $) = I{a^, $) + /(a(l - $). 

We now show that if a G S^, for some r, then /(a(l — $) = 0{k^^) for any i? > 0. The 
implicit constant is of linear growth in r. 

In the usual way, we integrate by parts repeatedly with the operator: 



that is, we apply the transpose 



k NV$| 



to the amplitude. The second term is a scalar multiplication. 

We first observe that the coefficients of C belong to S'°(y x M x R). Indeed, by assumption 
L G S^(Y). In the expression 

, v<i> A<i) v2<i>(v<i>,v<i>) _ \L wlhiyiyL) 
' fv$|^ ~ |v$|2 ^ |v$|4 ~ |v$|2 ^ pv$p ' 

it is then obvious that the numerator belongs to S^{Y). Since the denominator |V$p = 
I VyLp is bounded below on supp(l — ip), it follows that the coefficients belong to S^{Y). 

We now verify that each partial integration lowers the symbol order by one unit of k~^ , 
i.e. that {k-^C')^A{k,y) G S^^\Y). We prove this by induction on i?. As i? — i? + 1, we 
apply one of two terms of L*. We know that each differentiation improves the symbol order 
by one unit of k"^ . But each coefficient multiplies by an element of 5'°(y), hence preserves 
symbol order. 

The remainder has the form 

(135) |i?emM,i?(a, < sup max - - Ai)^fc((t, /i, a;, 0)|. 

Since - - /i, a;, 0) G S°(M x R x F), the right side is 0{k^^~^'^^^). 

If we choose 5 satisfying 1 > 5 > 1/2, then Jz-R]^{^-S)2R _ j^-{25~-i)r gjygg negative power 
of k. This is sufficient for the proof of the remainder estimate. 
This completes the proof. 



□ 



Combining Lemmas 19 . 2119 . 31 we obtain 

- ?(02), (pljAikll, 01, ... , (t)M-\a\)p{t)dtdlld(t)i ■ ■ ■ d(j)M-\a\+l 

(136) 

g(ifc-riogfe)rL^ fc-l'^l{^f^^(a^.,,j + a^-.,,j)fc-^' + 0(A;-«)}, (M > mr), 
Onikr^), M <mr 
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Indeed, the critical point set of the phase $ = t(l — /i) + jj,C is given by: 
(137) 

= • • • ,0A/-|<7|) : /i = l,t = £(01, . . . ,0Af-a|),g(0l), • • • ,g(0A^-|a|) = l'']- 

Clearly, the stationary phase set is empty if M < mr. When M = rm the phase is non- 
degenerate and we obtain the result stated in Proposition 19.11 by stationary phase. This 
completes the proof of Theorem 11.11 
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